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problems related to the Fourier representation 
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value problems. A minimum number of results 
are proved and the student can concentrate 
exclusively on application of results without 
being side-tracked by lengthy proofs of 
subsidiary points. 

The solution of boundary value problems 
using series of orthogonal functions other 
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and the application of Fourier and Laplace 
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Chapter 1 
Fourier Series 


1.1 Piecewise Smooth Functions A function j{x) is said to be piecewise 
continuous in the interval [a, 6] (the notation [a, b] will be used to denote 
the closed interval a ^ x ^ b) if there exist a finite number of points 
x u x 2 ,...'>x n (a = x t < x 2 ... < x n = b) such that /(x) is continuous in 
each sub-interval x f < x < x f+1 and has finite one-sided limits at the 
end points of each interval; these left- and right-hand limits will be denoted 
by /(x.+0) and /(x i+l -0) respectively. Clearly / need not necessarily 
even be defined at the end points of the sub-intervals. 

A function / is said to be piecewise smooth in the interval [a. ft] if / 
and its first derivative /' are both piecewise continuous in this interval. 

Problem 1.1 Show that, in the following, /(x) is piecewise smooth in its 
interval of definition. Determine the points at which / or /' are dis¬ 
continuous. 

(a) /(x) = x+c, ( c ^ x ^ 0), /(x) = c x, (0 ^ x ^ c). 

(b) /(x) = 0, (-1 $ x ^ 0), /(x) = x, (0 ^ x ^ 1). 

Solution, (a) / is obviously continuous in each of the intervals 
— c ^ x < 0, 0 < x ^ c and /(0+0) = /(0 — 0) = c; thus / is con¬ 
tinuous for x in [-c, c]. For x < 0, /' = 1 and for x > 0, /' = — 1; 
hence/' is continuous in each of the sub-intervals — c<x<0,0<x^c, 
but is discontinuous at x = 0. /' is not defined at the origin but both 
/'(0+0) and /'(0-0) exist and are finite (being + 1 respectively). Thus 
/' is piecewise continuous and / is piecewise smooth in [—c,c]. The 
graph in Fig 1 represents /. 


f f 










(b) / is clearly continuous in each of the regions — 1 $ x < 0, 

0 < x ^ 1 and /(0 + 0), /(0-0) are both equal to zero. Also /' = 0 for 
— 1 ^ x < 0 and /' = 1 for 0 < x < 1, so that /' is discontinuous at 
x = 0. /' is not defined at x = 0 but both /'(0—0) and /'(0+0) exist and 
are finite, being equal to 0 and 1 respectively. Thus / is a piecewise smooth 
function. □ 

1.2 Fourier Series If a function f[x) is integrable over the interval 
[«,u + 2c] then the Fourier series of / on the interval is the trigonometric 
series 
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Fourier series are often defined over one of the intervals [0, 2n], [—7t,7t], 
the results for these are obtained from (1.2), (1.3) with appropriate choices 
of a and c. The a r and b r are known as the Fourier coefficients of /. 

[An alternative form for the Fourier series can be obtained by expressing 
cos rmx/c, sin nnx/c, in terms of exp (± innx/c) and this gives 

OO 00 

a 0 + X a (i exp(in7tx/c)+ X a„exp(— innx/c), 

n = 1 n - 1 

where the bar denotes the complex conjugate and 

J Ca + 2c 

= 2c J ^ CXP * “ innx ^ dx ' 

This form can be further simplified by replacing n by —n in the second 

oc 

series, giving £ a n exp(hmx/c)]. 

n - — oo 

Theorem 1 If f(x) is piecewise smooth in the interval [n,n + 2c] then 
the Fourier series of / on [n,a + 2c] converges for all x in the interval. 
The sum, S(x), is equal to /(x) at all points where / is continuous; if / 
has a discontinuity at x = x r then 2S(x r ) = /(x r + 0)+/(x r -0). At the 
end points 2 S takes the value /(a+ 0)4* /(a 4-2c — 0). 


Theorem 2 (ParsevaTs 

[a, a 4* 2c] then 


theorem ). If / is piecewise continuous in 
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f 2 dx=\a 2 0 + X ia 2 +b 2 n ). 

n= t 


Problem 1.2 (i) Find the Fourier series of /(x) = x on [0,2zt] and 
determine the range of values of x for which the Fourier series converges 
to /(x). 

(ii) Sketch the function represented by the Fourier series in [-27t,47t]. 

Solution, (i) The coefficients a n and b n are found by setting a = 0 and 
c = n in (1.2) and (1.3). 

*2* j |*2ii 

x cos nx dx, b = - I x sin nx dx, 

71 Jo 


‘■41 

‘• 4 i: 


x dx = 27t. 

The integrals can be evaluated by integration by parts, 
f x 2 * 1 f 2n 

a n = I — sin nx I-sin nx dx, n # 0, 

|_7ro ^ = 0 nn Jo 

= -4— sin 2nn = 0, 
n 1 n 

[ x “i* = 2* i r 2 * 

— cos nx H-cos nx dx, 

nn Jx = o n7t Jo 


2 1 2 

= -h -j- sin 2nn = - - . 

n n z n n 


Hence setting c = it in (1.1) shows that the Fourier series of x in [0. 2n] is 



sinnx 
n * 


(ii) x is continuous and differentiable throughout the interval [0,2 tt] 
and hence, by Theorem 1, the Fourier series converges to x at all interior 
points. The value of the sum at the end points of the interval is the mean of 
the limiting values of x, and so the sum at each end point is n. The Fourier 
series has period 2n [i.e. S(x) = S(x + 27i)] because each term has this 
period. Thus the sum S(x) represents a function of period 2n, equal to 
x in 0 < x < 2n and to n at x = 0, 2n. S is shown in Fig. 2, where dots 
denote values at x = — 2n, 0, 27c, 47r. □ 


Generally, if S(x) = f(x) for 0 < x < 2c, then, outside the interval, 
S(x) will represent that function of period 2c which is equal to J (x) on 
0 < x < 2c and the Fourier series can be regarded as representing. 
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outside the interval 0 < x < 2c, the periodic extension of /. The Fourier 
series cannot represent a function f(x) for all x unless / has period 2c. 


Problem 13 Find the Fourier series of cos 4nx/c in [0,2c]. 


Solution. In this case a in (1.2) and (1.3) is zero and hence 

1 C 2c mix 4tcx . 
a = - cos - cos - ax 

" c Jo c c 

1 f 2c f (n - 4 )nx (n 4* 4)nxl 

-fcl L cos —— +0 “—r-J*- 

, if 2 *, mix 4nx 

b = - sin - cos - ax 

" c J 0 c c 

1 f 2c f . (n-4)nx . (n+4)jtx"| 

= — I sin--— + sin- - —\dx. 

2c <0 L c c J 

For n # 4 the integrals are all zero and hence a n = b n = 0, n # 4. Setting 
n = 4 gives a 4 = 1 and b 4 = 0, i.e. the Fourier series of J{x) = cos 4 nx/c 
is cos 4nx/c. □ 


This result illustrates an important theorem: If two trigonometric series 


/ mix , . nnx\ V / , nnx ,, . nnx\ 

l[ a n COS — + b n Sm —b L { a n COS — + b » S ' n —)> 

n- 0^ ' n = 0 ' ' 


converge to the same sum at all but at most a finite number of points in 
[a,a + 2c] then a n = a' n , b n = b' n . An equivalent statement is that if a 
Fourier series is identically zero then so are all the coefficients. 

Another corollary is that the Fourier series of any finite series of the form 

N 

Z l a n cos (nnx/c) 4- b n sin (nnx/c)] is the series itself. 

n = 0 
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This latter result can be proved by direct calculation of the Fourier 
coefficients, the coefficient of cos mnx/c in the Fourier series is, by defini¬ 
tion. 


1 

c 



a cos 


mix , . nnxl 

- + b„ sm- 

c " c 


mnx _ 

cos- dx. 

c 


The following results can be verified to be valid: 



nnx mnx 

cos -cos- 

c c 


dx 



. nnx mnx 

sin-cos- 

c c 

. nnx . mnx 

sin -sm- 

c c 


dx 


dx 


0, n # m. 



nnx 


cos 2 - dx 


c, n t* 0, 
2c. n = 0. 



sm 


, nnx , 
z — -dx 


jc, n t* 0, 

jo, k = 0. 


Since the summation is over a finite number of terms the order of sum¬ 
mation and integration can be interchanged and the above expressions 
show that the only non-zero contribution is that for m = n giving that 
the Fourier coefficient of cos mnx/c is a m ( m ^ N) and zero otherwise. 

A direct consequence of the above analysis is that if /(x) is assumed to 

X 

have an expansion of the form £ (a n cos(nnx/c)+b„sin(nnx/c)) then 

n = 0 

carrying out the above analysis (formally replacing N by oo) shows that 
a n and b n are the Fourier coefficients. This analysis is only formal as it has 
assumed a particular form for /(x) and also assumed that integration 
and summation can be interchanged in an infinite series. 

Problem 1.4 Find the Fourier series on [—1,1] of the function fix): 
f(x) = 0,-1 ^ x ^ 0; /(x) = x, 0 ^ x ^ 1. Hence show that 

f (2r-l)- 2 => 2 . 

r = 1 

Solution. In this example a = — landc = 1. The function /(x) does not 
have the same form over the whole range in the integrals defining a n and 
b n . Thus the interval of integration has to be split into the two intervals 
[-1,0] and [0,1]: 

a n = j* 1 ^ f{x)cos nnx dx = j*' x cos H7tx dx, [/= 0,-1 ^ x ^ 0] 

b n = J ' ^ f(x)smnnx dx = xsinM7ix dx. 

For n # 0 the integrals can be evaluated by integration by parts; the 
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case n = 0 involves an elementary integral giving a 0 = Thus 

fx sin nnx~\ x “ 1 if 1 . [coswr —1] [(— 1)" — 1 ] 

a = --sin nnx dx = -—j“ 2 —- = - j-j — 

L nn J* = o nn Jo n 71 n n 


|"x cos mtx~\ x “ 1 1 f 1 

L nn l-o nn Jo 


cos nnx dx 


nn 
cos nn 


nn 


-•4- — 2 — 2 s^n nn = 
n n 


(-1)' 


,n+ 1 


nn 


Thus a n = 0 (n even), and a 2r _ 
Fourier series is 


, = -2/(2r - l) 2 * 2 , r = 1, 2, 


The 


f 

4 7t 2 L 


cos(2r— l)7cx 


JU 

z 


( —l) n+1 sin nnx 


(2r — l) 2 nn 

Since / is a piecewise smooth function [Problem 1.1(b)] it follows from 
Theorem 1 that the sum of the Fourier series at x = 1 is 
2 [/(l)+/(—1)] = 2 * Substituting x = 1 in the Fourier series gives 

I (2r-l)- 2 = £* 2 . □ 

r = 1 

Problem 1.5 Find a Fourier series converging to |(sin x)| for every x. 

Solution. A Fourier series cannot represent a function for all x unless 
it has the same period as the function (in this case n\ see Problem 1.2). 
The precise interval is immaterial as both |(sinx)| and the relevant sines 
and cosines have the same period, and the integral of a periodic function 
over a cycle is independent of the end points. This can be proved as 
follows: Given F(x) = F(x + 2c) then (1-2), (1-3). Hence 

j fl + 2c F(x)dx = j 2c F(x)cix+ 2c F(x)dx = 

= f 2c F(x) dx + [.* +2c F(x + 2c) dx 

J a J 2c 

= j # 2c F(x) dx+ F(x) dx = J 2f Fix) dx. 

Thus the integral is independent of a. In the present problem we choose 
the interval [0, 71 ] for simplicity, i.e. a = 0, c = j 7 i in (1.2), (1.3). Hence 


-if. 

. 2f« . 

0 = - si 

" *J 


sin x cos 2nx dx = — 


7i(4n 2 -l) 


sin x sin 2 nx dx = 0. 

By Theorem 1 the Fourier series of |sinx| converges to (sinx| at all 
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interior points of [0,7r]. At the end points the series converges to the 
mean of the end values of | sin x|, i.e. to zero. It then follows by the periodi¬ 
city of | sin x | that, for all x. 


[sin x| 


2 4 cos2nx 

n n Li 4n 2 — 1 * 

n = 1 


□ 


Problem 1.6 If /(x) = sinx, 0 ^ x $ 7r, f(x) = 0, n ^ x ^ 2n, find the 
Fourier series of /(x) on [0,27r]. 


Solution. This problem can be solved most easily by noticing that 
/(x) = 3 (sinx + |sinx|). The Fourier series now follows directly from 
Problem 1.5 and the fact that sinx is its own Fourier series. The required 
Fourier series is 


^sinx-f 


i 


cos 2 nx 
(4 n 2 — 1)’ 


□ 


Problem 1.7 Find a Fourier series converging to sin 2 jx for all x. 


Solution. sin 2 jx is a function of period 27r and, as in Problem 1.5, it 
suffices to determine its Fourier representation in [0,27r]. This representa¬ 
tion can be obtained most simply by writing sin 2 ^x as ^(1 — cosx) and 
noticing that, by the theorem quoted in Problem 1.3, the Fourier series 
of this function is itself. □ 


1.3 Fourier Cosine and Sine Series In certain circumstances the Fourier 
series of a function /(x) in [ — c,c] consists only of sine terms or only of 
cosine terms. This arises when 

(i) /(x) is an even function of x. [i.e. /(x) = /( — x)]— cosine series 

(ii) /(x) is an odd function of x. [i.e. f(x) = — /( — x)]— sine series. 

This is illustrated in the following problem. 

Problem 1.8 (i) Show that the Fourier series on [ — c,c] of an even 
function /(x) contains no sine terms, (ii) Determine the general form of 
the Fourier series when /(x) is an odd function. 

Solution, (i) Setting a = —c in (1.3) gives 

cb n = J /(x) sin dx. 

The range of integration is now split into two ranges [ — c,0] and [0,c], 
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and x is replaced by — x in the former. Thus 


C K 


f° v . MIX _ f C . . MIX . 

f(x) sin - </x + /(x)sm - dx 

J - c c Jo c 

f C N . MIX , f C v . W7TX , 

— I /( — x)sin - dx+ /(x)sm - dx = 0, 

Jo c Jo c 


since /(x) = /( — x). Therefore the Fourier series will not contain any 
sine terms. From (1.1) we note also that 


ca n 


f ° v mix , f c x mix , 

/(x)cos - ax + /(x)cos - dx 

J -c C Jo C 

f C H7TX f C MIX C c 

/(-x)cos - dx + /(x)cos - rfx = 2 

Jo C Jo C Jo 


/(x) cos 


MIX 

c 


dx. 


(ii) For the case when / is odd, it follows from the above expressions 
for a n and b n , by replacing /( —x) by — /(x) that a n = 0, 

cb n = 2 Jq /(x) sin (nnx/c) dx , whence the result follows. □ 


Problem 1.9 Show that, in [ — 7r, xc]. 


X 2 = j7T 2 + 4 ") 


(— 1)" cos nx 


n = 1 


and hence evaluate (i) X n 2 ,(ii) X (— l) n+l n 2 , (iii) X 71 4 

n = 1 n = 1 n = 1 


Solution. The first step is to determine the Fourier series of x 2 on 
[ — 7i, 7r]. Since x 2 is even the series will consist only of cosine terms, and 
setting c = n in the appropriate equation of Problem 1.8 shows that the 
Fourier series is 


oo 2 f n 

? a o + X a n cos nx where a n = - I x 2 cos nx dx. 

«= i n Jo 

The integral can be evaluated by integration by parts for n # 0, giving 
n 2 a n = 4(- If; the integration for n = 0 is trivial and gives 3a 0 = 2k 2 . 
Since x 2 satisfies all the conditions of Theorem 1 and takes the same value 
at the end points of the interval [ — 7 r , 7 r ], the sum of the Fourier series is 
equal to x 2 at all points in [ — 7 t , 7 r ], i.e. 

x 2 = ^n 2 4-4 X [( — l)"cosnx/n 2 ], |x| ^ n. 

n =* 1 

T 

(i) Setting x = n gives, since cos wr = (— 1)", £ n~ 2 = jn 2 . 
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(ii) Similarly setting x = Ogives £ (-l)" +, n 2 = -fen 2 . 

n = 1 

oo 

(iii) The determination of £ n~* is less straightforward, but can be 

n as 1 

made by noticing that the terms are squares of the coefficients of the 
Fourier series for x 2 . This suggests use of Theorem 2 (Parsevafs theorem) 
which gives 



Hence 90 X = * 4 □ 

n = 1 

For any function /(x) defined in [0,c] one can construct an even 
function h(x) equal to f(x) for 0 < x < c by setting h(x) = /(-x) for 
— c < x < 0 and thus, by Problem 1.8, there will exist a Fourier series 
representation of /(x) in 0 < x < c consisting only of cosine terms. 
Similarly by defining h{x) by — /( — x) for — c < x < 0 it will be possible 
to construct a Fourier series representation of/(x) in 0 < x < c consisting 
only of sines. These representations are formally defined as follows: 


Fourier cosine series If/(x) is integrable on T0,c] then the Fourier cosine 
series of / on [0, c] is defined to be 


where 


2 - 0 + 2 , ‘ , . cos — 

n - 1 

2 f c „ , nnx 
a „ = -j o /(x) c ° s -d. 


(1.4) 


Fourier sine series lf^(x) is integrable on [0, c] then the Fourier sine 
series of / on [0, c] is defined to be 

, . nnx 

b n sin-, 

" c 

n - 1 

. 2 f c v . nnx . 

where b n = - /(x)sm- dx. (1.5) 

c Jo c 

Theorem 3 The Fourier cosine and sine series of a function /(x) which 
is piecewise smooth in [0, c] converge to /(x) at any point x in 0 < x < c 
where / is continuous. At any interior point x = x r , where / is discontin¬ 
uous, both series converge to \[f(x r + 0)+/(x r -0)]. The cosine series 
converges to /(04*0) at x = 0, to /(c — 0) at x = c and is always con¬ 
tinuous. The sine series converges to zero for x = 0, c. Theorem 3 can be 
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proved by on“ 

mSelval whose length is half the period of the functions involved. 

is that for a funcUon 

tSdwTtf* additional conditions /(0) - /(c) - 0. Another corollary 
function continuous in [0,c] and vanish,ng at* = Oand * - c 
can be represented throughout the interval by its Fourier sine senes. 

P,M,m 1 10 Find, on the interval [0,n], (a) the Fourier cosine seriB 
and (b) the Fourier sine series of cost*. Sketch the functtons defined 
by the Fourier series for |x| sS 3*. 

Solution, (a) The Fourier cosine series of x is \a 0 + ^ a „ cos " x w,hete 

= I2ln) f n x cos nx dx. We find a 0 = 7 t, for n # 0 integration by parts 

? 4 ° off-ii"-n In Thus a will vanish for even values of n. For 

gives n a n = 2|_(-1) U/*- ", n 2 a =_ 4. The Fourier 

„ = 2 r— 1 (r = 1,2,...) we get it(2r-l) a 2r _ l 

cosine series is thus 

t 4 £ cos(2r — l)x 
Zi (2r— l) 2 

r = 1 

and bv Theorem 3 converges to x in [0,7t]. , . rn . 

The cosine series is an even function of x and is equa to x m [ »*0 
- x in [ - n, 0], This even continuation (denoted by heavy lines) of x is 
shown in Fig. 3 together withits periodic extension (period 2ji)to [ - 3*. 3*]. 


(b) The coefficient b n in the sine series of cos^x is 
2 f n 1 f n 

b n = - sinnxcos^x dx = - [sin(n+^)x + sin(n —i)x] dx 

71 J o 71 J o 

8 n 

7r(4n 2 — 1)’ 

The series converges to cos^x at all interior points of [0,7i] and to zero 
at the end points and hence represents cos^x in 0 < x ^ n. 

The series defines an odd function of x and hence in — n ^ x < 0 is 
equal to — cos^x. From this and the periodicity property it follows that 
the graph is as shown in Fig. 4 (the form of the series in [-7i,7r] being 
denoted by heavy lines). □ 



Fig. 4 


Problem 1.11 (i) Determine the general form of the Fourier cosine and 

sine series in [0,c] of a function /(x) satisfying the relationship f(c — x) = 
/(*)• 

(ii) Obtain the corresponding series when /(c—x) = —/(x). 

Solution, (i) We consider first the integral from x = 0 to x = c of an 
arbitrary function F(x) such that F(c-x) = F(x). Such a function is 
said to be even about \c as setting x = \c—e gives F(^c+e) = F(^c—r.), 
i.e. it takes the same value at a given distance on either side of x = ^c. 
[If F(c —x) = — F(x) then the function is said to be odd about x = ^c.] 






















The range of integration is split into the two intervals [0, \c], [^c, c], thus 

f F(x)dx = \ iC F(x) dx+ 1° F(x)dx 
J o Jo J*c 

= j* c F(x)dx + j* c F(c —x)dx = 2 F{x)dx. 

For a function G(x) odd about x = \c it follows that G(x) dx = 0. 
Sin nnx/c is odd or even about \c according as to whether n is even or 
odd and cos nnx/c is odd or even about \c for according as to whether n 
is odd or even. Thus, for the given function f(x\ /(x)cos nnx/c is even or 
odd about \c according as to whether n is even or odd and from the above 
general result it follows that a n = 0 in the cosine series for n odd and 

4 r* 2nrx 

a 2r = - \ /tocos - dx , r = 0,1,.... 

c Jo c 

f(x) sin nnx/c is odd or even about \c according as to whether n is even or 
odd and thus b= 0 for n even and 

n 

4 P* c nx 

b 2 r-i = “ /(x) sin (2r — 1)— — dx, r = 1,2,.... 

c J o c 

(ii) In this case /(x) cos nnx/c is even or odd according as to whether 
n is odd or even and thus a n = 0 for n even and 
4 f± c jix 

a 2 r-i = - /(x) cos (2r — 1) — dx, r = 1,2,.... 

C J 0 c 

Similarly b n = 0 for n odd and 

4 [*** 2rnx 

b 2r = - fix) sin-dx, r = 1,2,.... □ 

c J o c 


Problem L12 Find the sine series for cosx in [0,7r]. 

Solution. Cosx is odd about \n and hence, from Problem 1.11, the sine 

00 

series will be of the form £ b 2r sin 2rx where 


■w-sjr 


cos x sin 2 rx dx = 


8r 


Jt(4r 2 -1) 


The sine series converges to zero at x = 0 and x = 7 t and hence 


cos x 




r sin 2rx 


0 < x < n. 
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Replacing x by ^x gives 


cos ;X 


5 y ^sinrx 
n L 4r 2 —T 


0 < x < 2n. 


This latter result is in agreement with one previously obtained in Problem 

1 . 10 . 

Problem 1.13 Find the sine and cosine series on [0,c] of the function 

f(x): fix) = x, 0 ^ x ^ \c: fix) = c —x, \c < x ^ c,and hence evaluate 

00 

Y, (2r— l) -2 . Sketch the functions defined by the Fourier series in 


\M 


^ 2c. 


Solution. The cosine series is \a Q 4* X a n cos i nnx /c) where 


n = 1 


2 f c _ / . nnx , 

a n = cJ 0 /( X)COS— ^ 


fix) does not have the same form over the whole interval of integration 
and this interval has to be split into the two intervals [0,^c], [yc,c]. / is 
even about x = \c and from Problem 1.11 it follows that a n = 0 for odd 
values of n. However in order to illustrate the point that the correct 
Fourier series can be calculated without direct use of the evenness property 
the integral for a n will be evaluated directly. 

ca f* nnx _ f c , , nnx . 

— 5 = xcos- dx 4- (c — x)cos-dx, n ^ 0 

2 jo C J*c C 

= j^2cosy-l-(-irj, n / 0 (by integration by parts). 


Clearly a n vanishes for n odd and hence 

a 2 r = (cos nr-1). 

This vanishes for r even and the only non-zero values are for r = 2s — 1 
[s = 1,2,...] when a As _ 2 = —2c/n 2 (2s— l) 2 , the Fourier series is thus 

1 2c y cos(4s—2)7tx/c 

4 ° it 2 4", (2 s- !) 2 

/(x) is piecewise smooth throughout [0, c]. Hence, by Theorem 3, the 
Fourier series converges to / on [0, c] and setting x = 0 gives 

£ (2s-l) -2 = |jt 2 . 

s= 1 
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The cosine series represents an even function of x equal to f(x) in [0, c] 
and its form is shown in Fig. 5. For c = 2n this has the same form as the 
Fourier series of Problem 1.10 (a). 

The coefficients b n can also be calculated directly as above or by using 
the evenness property of f(x). Problem 1.11 shows that b n = 0 (n even) and 

r ^ 7ix 

cb 2r _ l = 4 xsin(2r-l) — dx, r = 1,2. 

Jo c 

4(_ iy -1 c 2 

= —tttT’ integration by parts]. 

(2i— l)*ir 

f(x) vanishes at x = 0 and x = c and thus, in [0,c]. 


n 2 U(7r-\y c 


r = i (2r — 1) 


The sine series represents an odd function of x equal to /(x) on [0, c] and 
of period 2c: its form is shown in Fig. 6. 
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1.4 Cosine and Sine Series of Period 4c In applications it is often 
necessary to represent a function /(x) on [0,c] as a cosine or sine series 
whose general terms are a 2n _, cos(2n — l)7tx/2c and b 2n _ x sin(2n — \)nx/2c 
respectively. From Problem 1.11 it is seen that such representations could 
be obtained on [0,2c] for functions odd or even about x = c, and thus all 
that it is necessary to do is to extend the given function to [0,2c] and make 
it even or odd about x = c as necessary. The coefficients a 2n _ x , b 2n _ x can 
be obtained by replacing c by 2c in Problem 1.11, giving 


-if 

-if 


(2n-l)nx 
J (x)cos-—-dx. 


f{x) sin 


2c 

(2m— 1)tcx 
2c 


dx. 


The even extension of / is continuous and thus the sine series converges 
to/inO < x ^ c and to zero at x = 0; the odd extension is discontinuous 
at x = c and the cosine series converges to / in 0 ^ x < c and to zero 
at x = c. 


Problem 1A4 Determine the Fourier cosine series and sine series of 
period 4n which represent x in [0,7i] and sketch the functions represented 
by the series in |x| ^ 4n. 

Solution. Setting c = n in the preceding analysis gives 

00 

X= X a 2n-l COS 2^2w- 1)X, 0 ^ X < 7t, 

n = 1 

x 

X= X ^2n-! sin 2(2n—1)** 0 < X ^ 7C, 

n = 1 

where 

2 f* 2 f n 

a 2n -\ = ~ xcos j(2m— l)x dx, b 2 . =- xsini(2n— l)x dx. 

71 Jo n Jo 

The integrals can be very easily evaluated by integration by parts giving 

4(-l)" +1 8 8( —iy ,+ l 

“ (2m — 1) Jt(2n — 1 ) 2 * 2n ~ i ~ n(2n-l) 2 ‘ 

x vanishes at x = 0 and hence the sine series converges to x in [0,7c]. 
Cos^(2m — l)(7r — x) = — cosj(2n— 1)(tc + x) and thus the cosine series 
represents a function odd about x = n (this of course also follows from the 
discussion preceding this example). Thus the form of the cosine series is 
known for 0 ^ x ^ 2n and, as the cosine series is an even function of 
x, it can be now determined for — 2n ^ x ^ 0. The series is thus known 
for a complete cycle and the form for all x can now be obtained and is 
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Fig. 7 


given in Fig. 7. The form for the sine series is obtained similarly and is 
given in Fig. 8. □ 



Fig. 8 


These results provide two more examples of Fourier series expansions 
which represent the same function over a given interval but which repre¬ 
sent very different functions when considered for all x (cf. Problem 1.13). 

Problem 1.15 Find the representation in [0,c] of sinnx/c as a Fourier 
cosine series of period 4c. 

Solution. The coefficient a 2n _ l of cos (2n—\)nx/c in the expansion is 
given by 

2 f c . 71X 1. 7tx , 

a 2n-i =~ ^sm-cos-(2n-\)—dx. 

The integral can be evaluated by expressing the integrand as a sum 
of sines and this gives a n = — 8/7i(2n + l)(2n — 3). □ 
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1.5 Differentiation and Integration of Fourier Series Differentiation of 
a Fourier series will produce another Fourier series whose coefficients are 
proportional to na h , nb n where a n and b n are the coefficients of the original 
series. Thus, even if the original series converges, the differentiated series 
may not even converge if a n and b n are not sufficiently small for large n 
and it is therefore useful to be able to derive, from the form of the original 
function, the behaviour of the coefficients for large n and we have the 
following theorems: 

Theorem 4 If / is piecewise smooth in [a,a + 2c] then a n and b n are at 
most 0(n~ l ) (i.e. lim na n , lim nb n are both finite). 

n~*co n-*oo 

Theorem 5 If f and its first (p— 1) derivatives are continuous in 
[a,a+2c]; f {r) (a) = f {r) (a + 2c\ r = 0, 1, (where / (r) is the rth 

derivative of /) and / (p) is piecewise smooth in [a,a + 2c] then a n and 
b n are at most O (n ~ p ~ l ). 

Problem 1.16 Show that, for large n, the coefficients of the cosine series 
on [0,7i] of a function /(x), continuous with a piecewise smooth first 
derivative, are 0(n“ 2 ). 

Solution. The behaviour of the coefficients cannot be deduced directly 
from the above theorems which refer to the complete Fourier expansion 
of a function. The cosine series expansion of / is the complete Fourier 
expansion on [ — 7r,7c] of the extended function even in x and equal to 
/ in [0,7r]. This function is continuous in [ — tt,7c], taking the same value 
at both end points, and hence Theorem 5 may be applied with p = 1, 
giving that the Fourier coefficients are 0(n" 2 ). □ 

The derivative of the extended function is not continuous at x = 0 as 
its left and right hand limits are + /'(0) respectively, at the end points this 
derivative is ±f'(n). Thus, if /'(0) = f\n) = 0, Theorem 5 shows that 
the Fourier coefficients of the cosine series for / would be 0(n~ 3 ) if /' 
were continuous and /" piecewise smooth. 

Problems 1.10(a), 1.13, 1.15 are all examples on deriving cosine series 
of functions satisfying the conditions of the present example and all have 
Fourier coefficients of 0(n~ 2 ). It follows that the coefficients of the sine 
series obtained by differentiating cosine series of the above function are 
0(n -1 ) and the series is at least convergent. This however does not prove 
that the differentiated series converges to /'(x). Similarly the coefficients 
of the Fourier sine series off are 0(n _1 ) but will be 0(n" 2 ) provided that 
/(0) = f(n) = 0. Thus in many cases the series obtained on differentiating 
the sine series of a continuous function will not even converge. 
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The general conditions concerning differentiation of a Fourier series 
are contained in the following theorem: 

Theorem 6 If /(x) is continuous and /'(x ) piecewise smooth in [0, c] then 

(i) The Fourier cosine series for / may be differentiated term by term 
and the resulting series converges to /' in 0 < x < c. 

(ii) The Fourier sine series may be differentiated term by term if, and 
only if, /( 0 ) = /(c) = 0 , and the differentiated series then converges 
to /' in [ 0 , c]. 

For the case when /(0) or f(c) are not zero the coefficients of the cosine 
series for /' are given by 

« 0 = \ a n = ^- 2 t [m+(-ir l f(c)i 

In applying this theorem to cases where the derivative of/ is not defined 
at a point x = x 0 , /' at that point is to be interpreted as 
2 [/ ( x o + ®) + / ( x o~^]- The theorem is proved by using integration by 
parts to express the Fourier coefficients of /' in terms of those for / 


Problem 1.17 Use the results of Problem 1.10(b) to obtain the Fourier 
cosine series of sin\x on [ 0 , ri]. 

Solution. Problem 1.10(b) shows that 


1 8 ^ 
cos- X = - ) 
2 n 


n sin nx 
4n 2 -1 ’ 


0 < x ^ n. 


The Fourier coefficients a n of — ^sin^x are, from Theorem 6 (ii), 


a o “ _» 


8 n 2 


a = 


n(4n 2 — 1) 


2 

n 


n(4n 2 - 1)’ 


_ . . 2 4 p cos nx 

Thus sin^x =-) — 5 —0 < x ^ n. □ 

z n njL* x (4n 2 — 1 ) 

Problem 1.18 Use the results of Problem 1.13 to determine the cosine and 
sine series in [ 0 ,c] of the function g(x): g(x) = 1 , 0 < x < ^c; g(^c) = 0 ; 
g(x) = -l,\c < x ^ c. 


Solution. Except at x = the function g is clearly the derivative of 
/ defined in Problem 1.13. At x = ^ c , however, /' is not defined and thus 
/' at x = for application of Theorem 6 , has to be interpreted as the 
mean of the limiting values of the derivative on either side of x = 
i.e. zero. Theorem 6 (i) then shows that the sine series for g is 

ao 

4]T [sin(4s - 2)nx/c]/{2s - 1). Theorem 6 (ii) shows that the coefficients 
1 
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in the cosine series of /'are given bya 0 = 0[/(0) = /(c) = 0], a n - nnbjc , 
where 4(—l) r ” l c 

b n = 0 (/i even), b 2 r _i= {2 r-\) 2 i? % r = 

This gives a n = 0 (n even) and a 2r _ x = 4(- l) r_ 1 /(2r— l)7r,r = 1,2.... □ 

Results analogous to those of Theorem 6 can be obtained for the cosine 
and sine series of period 4c defined in section 1.3. 

Theorem 7 If / is continuous in [0,c] with a continuous first derivative 
then 

(i) when / is represented in [0, c] by the series 

* 

TtX 


Yj fl 2,-l C “( 2r -% 


r = 1 

the derivative /' can be represented by the series 

oc 

nx 


where 


Y b 2r— iSin^r-l)^ 

f = 1 

2(-l) r ' M 


2f — 1 


V (2r—1) 

f(c)-~^na 2 r _ r 


(ii) when / is represented in [ 0 , c] by the series 


r Ib 2r _ 1 sin(2r —1)^ 


the derivative / # can be represented by the series 

OO 

Y Vr' cos ( 2r - 1 )| 


where 


r = 1 


*2,-1 = -“/(PH 


( 2 r-l) 7 tb 2f _ 1 

2 c 


This theorem is also proved by means of integration by parts. 


Integration of Fourier series The situation concerning the integration 
of Fourier series is much simpler than that for differentiation and it can 
be proved that a Fourier series of a function / may be integrated term by 
term and that the integrated series converges to the integral of f. 

Problem 1.19 Show that the Fourier cosine series of sin \x on [0, rc] may 
be obtained by integration from Problem 1.10. 


Solution. In Problem 1.10 it is shown that 


cos^f = 



n 

T~i —T sm nt ' 
4n — 1 


0 < t ^ n. 
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and integrating this identity from t = 0 to t = x gives 

... 8 (cosnx—1) 

2sm '* - L 4,‘-i ■ 0 «* s ”• 

n = 1 
oc 

It is known that £ (4n 2 —1)~ 1 = | and thus 


n- 1 


. , 2 4 ^ COS MX 

sin^x =-> —^— . 

2 n n Lt 4n 2 —1 

n = 1 


r 

The sum £ (4n 2 — 1) 1 can be calculated by setting x = 0 in the co- 
L n= 1 


sine series for sin x in [0, n \: 


2 4 ri cos 2 nx 

sin x =-> —--. 

n n Lt 4n 2 — 1 

n= 1 


0 ^ X ^ 


-] 


For a 0 # 0 the integrated series will not be a Fourier series. 


EXERCISES 


1. Show that 


r\ 2 v< cos/jx v-> sinnx 

(,)X = — + 4 Z^- 47t S— ^ < x < 2n. 


n = 1 


(ii) expax = 2 sinh H l + f (-iracosnx-nsinnx l 
n I 2 a a 2 + n 2 

*“ n- 1 - 1 

— 71 < X < 71. 

By considering the value of the series at the end points show that 

j ® 2 

2 na coth Tta = - + y , , a . . 

2 2 L, (a 2 + n 2 ) 


n = 1 


(iii) ^2-x(7i—x)(27t —x) = f 0 «; x s? 2 n. 

Lu n 
n = 1 


/• v i | . 4 ^ cos( 2 n-l)x 

00 00 

and hence evaluate £ (2n-l) -2 , £ (2n-l)~ 4 . 


n = 1 


n = 1 


2 . Determine on the interval [0,7t] the sine series for x and the ‘cosine’ 
series for cos ^x and obtain the range of values of x for which the series 
converge to the relevant functions. 
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3. Obtain by direct calculation the Fourier cosine and sine series in 
[0,c] of the function g(x) defined in Problem 1.18. 

4. Find the cosine and sine series of period 4c representing expx in 
0 < x < c. 

5. Determine the cosine series in [0,7t] for sin x by applying Theorem 6 
to the sine series for cos x derived in Problem 1.12. 
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Chapter 2 

Solution of Boundary-Value Problems by Means 
of Fourier Series 

2.1 Introduction Most of the subsequent problems will be concerned 
with the solution of the simpler partial differential equations of mathemat¬ 
ical physics. Particular examples of these equations are: 

V 2 w = p (Poisson's equation), V 2 u = 0 (Laplace's equation), 

- 1 d 2 u . . v . . 

v u = ~pJT (wave equation), V*u = k— (diffusion equation), 

V 2 u = ku (Helmholtz’ equation). 

For two independent variables (which is the case with which we shall 
mainly be concerned) all the above are particular cases of 

+ a 3 + b i Oofp!++ b 3 (y)u = f(x,y), ( 2 . 1 ) 


where a t , b t are known functions of x and y respectively and / is another 
known function of these variables. 

The methods to be described are only suitable for solving (2.1) in a 
rectangular domain D of the x, y plane and in the present chapter it is 
necessary that the range of one independent variable is finite and that the 
coefficient of the second derivative with respect to that variable does not 
vanish within D. Without loss of generality x is taken to be that variable 
(i.e. a, ^ 0 in D) and D to be the domain 0 x ^ c, b < y < d; either 
or both of b, d can be unbounded. 

The boundary conditions under which (2.1) is soluble depend critically 
on the sign of b l /a l (it will be assumed this sign does not change within D). 
[An equivalent method of relating the boundary conditions to the 
equation is by reference to the classification of partial differential equations 
in two independent variables (cf. P. M. Morse and H. Feshbach, Methods 
of Theoretical Physics. McGraw-Hill. 1953, Ch. 6 ). In this nomenclature, 
b l /a l being positive, negative or zero determines whether ( 2 . 1 ) is elliptic, 
hyperbolic or parabolic.] 

The boundary value problems to be solved in subsequent problems 
will be such that they possess unique solutions and the following table will, 
for most practical situations, indicate the appropriate conditions to be 
applied. 
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bja v + ve 
bja v -ve 
bja v 0 


u or cu/cy given on y = b.y = d. 
u, du/dy given on y = b. 
u given on y = b. 


u or cu/dx given 
for x = 0 . c. 


2.2 Homogeneous Problems Homogeneous problems in the present 
context are defined to be those of solving 


V* 




8 2 u 


du 


dx 2+ ° 2 dx 


d 2 u du 

a 3 U + h l-^2+ b 2fy + b 3 U = 0 . 


( 2 . 2 ) 


in D with u or du/dx vanishing on x = 0, c. The use of Fourier series (and 
their generalizations) to solve (2.2) relies on two separate concepts. The 
first is the principle of superposition which states that if any number of 
functions separately satisfy ( 2 . 2 ) then so will any linear combination ofthese 
functions (this is because the left hand side is linear in u). The second con¬ 
cept is that of separation of variables which consists of seeking product 
solutions of (2.2) of the form X(x)F(y). The basic technique in subsequent 
problems is essentially the determination of combinations of product 
solutions which are such that they satisfy all the conditions imposed. 
This approach is generally known as Fourier’s method and is illustrated 
in the succeeding problems. 


Problem 2.1 Solve the equation d 2 u/dx 2 + 6 2 u/8y 2 =0, 0 < x < c, 
0 < y < d. under the boundary conditions i/( 0 . v) = u(c,y) = 0,0 < y < d, 
u(x. 0) = J(x), »(x, d) = 0, 0 < x < c, where (i) /(x) = sin 7 rx/c, (ii) 
/(x) = x(c-x). 

(One physical interpretation of this problem would be the determination 
of the steady state temperature in a rectangular region where three sides 
of the rectangle are kept at zero temperature whilst the fourth is main¬ 
tained at a known temperature /(x).) 


Solution. Assume u is of the form X (x) Y (y) where X and Y are functions 
of x and y respectively. Substitution of this in the equation for u gives 

XT _ Y" 

X ~ ~ T ' 

the dashes denoting derivatives with respect to the relevant variable. A 
function of x can only be equal to a function of y if both are constant. 
Hence X" = —XX, Y" = AX where A is a constant (the separation con¬ 
stant). The condition u(0,y) = u(c,y ) = 0 gives Af(0) = X(c) = 0 and the 
general solution for X is Ax + B, A = 0; Ccos X i x -I- D sin A*x, A / 0. For 
A = 0 the conditions X(0) = X{c) = 0 show that A = B = 0 and. for 
A ¥= O.A’(O) = 0 gives C = Oand X{c) - 0 then requires that D sin X*c = 0. 
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Thus in order to have non-trivial solutions (i.e. D ^ 0)cA i must be equal 
to nn . where n can be any integer. X will then be proportional to 
sin nnx/c and the equation for Y can be solved showing that Y will be an 
arbitrary linear combination of cosh nny/c, sinhn 7 ry/c. Thus possible 
solutions for u. satisfying the conditions at x = 0 and x = c, are 
[A n coshnny/c + B n s'mh nny/c^sin nnx/c where n is any integer and the 
arbitrary constants A n , B n may depend on n. 

In accordance with the principle of superposition we thus consider the 
general form 


«- i f 

n - i L 


A. cosh + B_ sinh — 


sin 


nnx 


and attempt to find A n and B n to satisfy the conditions at y = 0 and 
y = d. The condition on y = d will be satisfied if 


, nn( t „ . , nnd 

A n cosh -= - B m sinh- 

" c " c 


and, for case (i), the condition on y = 0 gives 


. tlx nnx 

sin — = ) A sin-, 

c Lj " c 

it - 1 


0 < x < c. 


Thus the A n are the coefficients of the Fourier sine series in [0, c] of sin nx/c; 
this latter series is the function itself [cf. Problem 1.3] and hence A x = 1, 
A n =*0, n ^ 1 and u = — cosech (nd/c) sinh n(y — d)/c sin nx/c. It can be 
verified directly that this satisfies all the conditions imposed. 

The only difference in case (ii) is that A n will now be the coefficients 
of the sine series representing x(c —x) in [ 0 ,c], i.e. 


A 


n 


2 C e nnx 4 c 2 

-I x(c—x)sin - ax = — s—s 11 — cos nn\ 

c Jo c n J n J L 


on integration by parts. Thus A n = 0 (n even) and A 2r _ t = [ 8 c 2 / 7 t 3 ( 2 r-1) 3 ], 
r = 1,2,.... The complete solution is thus 


8 c 2 y 1 
7 r 3 ( 2 r — l ) 3 




x sin 


(2r— l) 7 rx 


c 

The fact that each term of an infinite series satisfies a differential 
equation does not necessarily mean that the sum satisfies this equation, 
as this would imply that differentiation of the series produces another 
convergent series. This is not always the case [cf. remarks following 
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Problem 1.16]. In the above series, however, the coefficients tend to 
zero sufficiently quickly for the series to be differentiated twice and the 
sum does satisfy Laplace’s equation. The boundary conditions are all 
satisfied (by construction) and hence the above series represents the 
solution for case (ii). □ 

Problem 2.2 Solve the previous boundary value problem when 
f(x) = x(c-x) and the condition on X = 0 is replaced by u(0, y) = y(d-y ), 
0 < y < d. 

Solution. The boundary conditions on neither x = constant nor y = 
constant are homogeneous but, as Laplace’s equation is linear, u can 
be written as Uj+m 2 , where u x and u 2 both satisfy Laplace’s equation 
and satisfy homogeneous conditions in one or other of the variables. We 
have 

= «j(c,y) = 0 , 0 < y < d, 

u,(x, 0 ) = x(c-x), u,(x,<f) = 0 , 0 < x < c, 

u 2 (x, 0) = u 2 (x,d) = 0 , 0 < x < c, 

u 2 ( 0 , y) = y(d-y), u 2 {c, y) = 0 , 0 < y < d. 

u, is the solution in Problem 2.1, case (ii). The problem for u 2 is very 
similar to that for u, except that homogeneous boundary conditions are 
prescribed on y = 0 and y - d. A solution could be obtained by a very 
similar procedure to that of Problem 2.1 but in this particular case all 
additional work can be avoided by noticing that interchanging x and y 
and c and d reduces the problem for u 2 to that for u x . Thus the solution 
can be obtained from the previous result as a sum of two infinite series. 

This technique of splitting the problem into two separate ones with 
homogeneous conditions to be satisfied on one family of parallel lines 
in each case may be employed for any linear elliptic equation when 
non-homogeneous conditions are prescribed on the boundary of a 
rectangular region. □ 


Problem 2.3 Solve the equation 

d 2 u 1 d 2 u . A 

— 2 , 0 < x < /, t > 0 . 

dx 2 v 2 3t 2 

where v is a constant, under the conditions 

u( 0 , t) = 0 , du(x,t)/dx = 0 , x = /; t > 0 , 

u( 0 ,x) = x, du(x,t)/dt = 0 , t = 0 ; 0 < x < /. 

(This problem can be interpreted physically as determining the longi¬ 
tudinal vibrations of a rod of length / clamped at x = 0 and free at x = I. 
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Initially the longitudinal displacement at a distance x from the fixed 
end is equal to x and the rod is at rest. These initial conditions could be 
produced by applying a force at the free end parallel to the rod and remov¬ 
ing this force at time t = 0 .) 


Solution. Assuming a product solution X(x)T(t) and proceeding as in 
Problem 2.1 gives X " = —AX, T" = — kv 2 T. The conditions on X are 
X(0) = X'(t) = 0 and the general form of X will be that of Problem 2.1. 
For A = 0 the conditions X(0) = X\l ) = 0 again give A = B = 0; 
for A ^ 0, X(0) = 0 gives C = 0 whilst X'(t) = 0 requires that D cos A*/ = 0. 
Thus lor a non-trivial solution A*/ = (2n — l) 7 i/ 2 , n = 1,2,..., where n is 
any integer and the corresponding X is proportional to sin (2 n- \ )nx/2l. 
Hence solving the equation for T and using the principle of superposition 
suggests the general form 


-it* 

n = 1 L 


COS 


( 2 n- \ )nvt 


B sin 


(2n — l) 7 u;t~| . (2n—l) 7 rx 


n. i 

- sin- 


21 


21 21 

The condition du/dt = 0, t = 0, is satisfied by taking B n = 0 and 

00 

u(x,0) = x implies that x = £ -4 n sin(2n— l)nx/2l. Thus 

8/(— 1) B+I 




(2n —l) 7 tx . 
x sin -—- dx = 


and 


u — 


8 / 

_2 


V ^ ~ ^ 

Lj (2h — 


21 

,n+l 


( 2 n-l) 


cos 


n 2 {2n 
(2 n— 1 )nvt 


21 


sin 


l) 2 

( 2 n — l) 7 tx 
~2l ’ 


u satisfies the boundary conditions by construction, but as it is not 
permissible to differentiate twice under the summation sign though the 
separate terms satisfy the wave equation, it is not possible to show the 
series also satisfies the equation, u can be written as F(x — vt) + F(x + vt) 
where 



n = 1 


(-I )"* 1 . (2n-\)nx 

(2^1)2 Sin ~^r- 


F(x) is related to the function sketched in Fig. 8 and is in fact the sine 
series of period 41 representing x in [0, /]. By analogy with Fig. 8 it has 
a discontinuous derivative at x = ^/,|/,..., and its second derivative 
does not exist at these points. Thus, as any twice differentiable function 
of x±vt satisfies the wave equation, F(x±vt ) are solutions of the wave 
equation for a given t at all but a finite number of points. 

The solution obtained is a solution in a somewhat generalized sense 
and in the general nomenclature of hyperbolic equations is called a 
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‘weak' solution. The discontinuity in the derivative of F essentially 
reflects the discontinuity inherent in the initial conditions. □ 

Simplified procedure In both Problems 2.1 and 2.3 the equation for 
X reduced to the simple harmonic one and it can be seen from these 
problems that this equation will always be obtained when the variable 
x only occurs explicitly in the equation as a constant multiple of d 2 u/dx 2 . 
The same situation will also occur if there is a constant multiple of u in 
the equation. Also, by analogy with the previous problems, possible 
solutions X can be obtained for those cases when u or du/dx vanish at 
x = 0, c. The forms for two cases have already been found and the results 
for the other cases are 

(i) du/dx = 0, x = 0,c: X = cos nnx/c , n = 0,1,2,..., 

(ii) du/dx = 0 , x = 0 ; u = 0 , x = c: 

X = cos(2n- l) 7 rx/ 2 c, n = 1,2,.... 

Thus for equations of the form of (2.2) with a 2 = 0 and a v a 3 constant 
it is possible to predict, for certain boundary conditions, without going 
through the general separation of variables procedure, the general form 
of the solution. It is useful to have the appropriate representations in the 
form of a table. 


Table 1. Representations for case a 2 = 0, a t ,a 3 constant 


Boundary conditions 

Representation for u 

(i) u = 0 , x = 0 , x = c 

00 

£ V„(>)sin nnx/c 

n * 1 

(ii) du/dx = 0 , x = 0 , x = c 

oo 

£ y n (y)COS wttx/c 


n = 0 

(iii) u = 0 , x = 0 ; du/dx = 0 , x = c 

X ^(y)sin( 2 n-l) 7 tx/ 2 c* 

n = l 

(iv) du/dx = 0 , x = 0 ; u = 0 , x = c 

CO 

£ y„( v) cos (2 n — 1 )nx/2c 

1 


Entries (i) and (iii) respectively would have been appropriate for Problems 
2.1 and 2.3. Once the appropriate representation has been chosen from 
the table then substituting this in the relevant equation and equating the 
coefficients of the trigonometric functions to zero (cf. remarks in Problem 
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c 

























1.3) gives the equation for Y n . (Theorems 6 and 7 show that, provided u 
satisfies the relevant homogeneous conditions, the series can be dif¬ 
ferentiated term by term twice.) The approach is illustrated in the following 
problems. 

Problem 2.4 Solve the equation d 2 u/dx 2 = 8u/dt, 0 < x < I, t > 0, 
under the boundary conditions du(x,t)/8x = 0, x = 0; u(l,t) = 0 , t > 0; 
w(x,0) = sin(7rx//), 0 < x < l. (Physically this corresponds to calculating 
the temperature in the slab 0 < x < / when the end x = 0 is insulated 
whilst the end x = / is maintained at zero temperature, the initial tempera¬ 
ture being sin(7rx//).) 


Solution. Table l(iv) gives the appropriate representation as u = 


00 

£ T n (t) cos (2 n — l) 7 tx/ 2 /. Substituting this into the equation gives 

n = 1 


-l 


(2w — l) 2 7t 2 
4 1 2 


T n (t) cos 


( 2 n-l) 7 tx ^ dT n (2n-l)nx 

21 L dt C ° S 2 / ’ 


n = 1 n = 1 

equating the coefficients of cos( 2 n — \)nx/2l in this equation gives 
dT n (2n—l) 2 n 2 ^ 

dt ~ 4 l 2 

and the general solution is 7^ = /l n exp[ — (2n — l) 2 ;r 2 t/4/ 2 ]. 

The boundary condition at t = 0 shows that T n (0) (= AJ are the 
Fourier coefficients in the expansion of sin nx/l as a cosine series of 
period 4/. This was considered in Problem 1.15 and A n can be found from 
that problem by replacing c by /. Hence 


8 v 1 exp[-(2n —l) 2 7 t 2 r/ 4 / 2 ] ( 2 n-l) 7 tx 

u =-> -cos-. 

7 iL x (2n+l)(2n-3) 21 

The coefficients in this series are exponentially damped for large n. 
Thus the order of differentiation and summation can be interchanged and 
u satisfies the differential equation ;u also satisfies the boundary conditions 
(by construction). □ 


Problem 2.5 Obtain a solution u of Laplace’s equation finite everywhere 
in the sector O^r^a, 0 < 0 ^ a, where r ,0 are polar coordinates, 
and such that it has zero normal derivative on the lines 0 = 0 and 0 = a 
and is equal to 0 on r = a. 

Solution, u satisfies r 2 8 2 u/dr 2 +r du/dr+ d 2 u/dQ 2 = 0, which is of the 

form of equation (2.2) with x replaced by 0. Table l(ii) gives the appropriate 

00 

representation for u to be £ FJr) cos nnO/cL. Substituting this in the 

n = 0 
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equation and equating the coefficients of cos nnO/oi. to zero gives 
r 2 F'; + rF' n -(n 2 n 2 /x 2 )F n = 0. 

For n ^ 0, independent solutions of this equation are r ±nitl * and the 
finiteness condition shows that the positive sign must be chosen. For 
n = 0 the general solution is A\nr + B and for finiteness A = 0. The 
appropriate general solution for F n is thus A n r nn,a (all n) and the boundary 
condition on r = a shows that A n a nnl * are the Fourier coefficients of the 
cosine series of 0 in [ 0 , a], i.e. 

A 0 = - P e do = \oi, A n a nnl * = - P e cos — dd. 

0 “Jo 2 “Jo « 

Integration shows that A n = 0 (n even) and A lr _ l = — 4 a/ 7 r 2 ( 2 r—l) 2 , 
[r = 1,2, . . .] (cf. Problem 1.10(a)). The infinite series for the solution 
can now be written down. □ 

The above analysis may be generalized to the case where u satisfies 
the equation (V 2 + /c 2 )u = 0 where k is a real constant, and the boundary 
conditions are unchanged. The equation for F n is 

a 2 r 2 F"4*a 2 rF| | +(/c 2 a 2 r 2 — n 2 n 2 )F H — 0, 


or 


#1. 

dz 2 




= 0 , 


where z = kr. This is Bessel’s equation of order wr/a, and the solution 
finite for z = 0 is J nn/ Jz) which is the Bessel function of order nn/<x. Thus 
the solution of this second problem is found by replacing r "* ,a , a nn,a in 
the solution found earlier by J nn/ Jkr\ J nn Jka). 

Problem 2.6 Find a function satisfying Laplace’s equation, finite and 
continuous everywhere in the circle 0 ^ r ^ a, and equal to cos 3 0 on 
the circumference, r, 0 being polar coordinates. 


Solution. There are no boundary conditions on 0 = constant which 
would let us employ Table 1 but assuming a solution F(r)G(0) gives 



G has the general solution A cos A*0 + B sin A*0 and A can be found by 
noticing that G will not be continuous and single valued in [0, 27 t] unless 
A = n 2 (n any integer). The solution F n corresponding to A = n 2 can be 
obtained from Problem 2.5 on setting a = n\ the finiteness condition 
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shows that F n will be proportional to r". Thus the general form for the 
solution will be 

ao 

A 0 + £ r*(A cos n0+B sin nff). 

n = 1 

Setting r = a shows that A n , B n can be calculated in terms of the Fourier 
coefficients of cos 3 0 and the latter are most easily calculated from cos 30 = 
4cos 3 0-3cos0. Hence B n = A n = 0, n ^ 1,3, A x a = * A 3 a 3 = □ 

Problem 2.7 Find u satisfying V 2 u = du/dt , t > 0, finite everywhere 
within the sphere r < a, vanishing on r = a and equal to unity when 
t = 0 . 

(This problem is equivalent to the determination of the temperature 
in a spherical region, the boundary of which is maintained at zero 
temperature, the initial temperature being prescribed.) 

Solution. The boundary conditions are spherically symmetric and it 
is reasonable to assume that u is a function of r only, hence 

1 d ( 2 cu\ du 

?s( r - F 

Writing ru as v gives d 2 v/dr 2 = dv/du and in order that u be finite when 
r = 0 it is necessary for v to vanish at r = 0 . v will also vanish on r = a 
and. as v satisfies an equation where r only occurs in the term d 2 v/dr . 

Table l(i) shows that the appropriate representation for v is 
00 

£ T n (t)sin rntr/a, and substituting this in the equation gives 

n = 1 

T' n = ( n 2 n 2 /a 2 )T n , with a general solution T n = A n t\\>(-n 2 n 2 t/a 2 ). The 
condition at t = 0 shows that A n are the Fourier coefficients of the sine 
series of r on [0, a], i.e. A n = (2/a) r sin(« 7 tr/a) dr = 2( - l f +l a/nn. □ 

Problem 2.8 Solve the equation d 2 u/dx 2 + d 2 u/dy 2 = d 2 u/dt 2 , t > 0, 
0 < x < c, 0 < y < d, under the conditions u = 0, x = 0, x = c; 
y = 0,y = d: u = xy(c-x)(d-y), du/dt = 0, t = 0,0 < x < c,0 < y < d. 

(A physical interpretation of this is the determination of the shape of 
a rectangular membrane given its initial shape and that its initial velocit> 
is zero.) 

Solution. As this problem involves three independent variables Table 1 
is not directly applicable, but assuming a product solution gives 

y» i j*rt 

~x + Y = T' 


30 


Both sides of this equation must be constant i.e. 



A. 


As X and Y are functions of independent variables x and y , both the 
terms in the first of the latter equations must be separately constant, i.e. 


XI 

Y 


= 



A. 


The equation for Y is similar to the equation of Problem 2.1 and in 
order that u vanishes at y = 0 , d it is necessary that n = m 2 n 2 /d 2 (m an 
integer). Similarly from the equation for X and the conditions at x = 0 . c 
it is found that (A -fi) = n 2 n 2 /c 2 (n an integer). Substituting A, fi in the 
equation for T shows that the general form for u is 

2, cos t+B"" sin f] sin — sin —, 

n = lm= 1 

where K 2 n = n 2 n 2 /c 2 + m 2 n 2 /d 2 . du/dt will vanish for t = 0 if B mn = 0 
and the other condition at f = 0 gives 

, vm * . . tmx . mny 

xy(c-x)(d-y) = ^ ^ ^ sin ~ sin ~y» 0 < x < c, 0 < y < 

() = lm = 1 

This is a double Fourier series and, for a general left-hand side. A mn 
could be found by successive application of equation (1.5), i.e. 

Amn = cd\\ xy(c ~ x ^ d ~ y)sin dxd y- 

This elaborate approach can be avoided in this case as the integrand is a 
product of a function of x and a function of y. The sine series for x(c — x) 
has been obtained in Problem 2.1 and that for v(d-y) can be derived 
from this by replacing x,c by y,d. Hence A mn = 0 if either m and/or n 
is even and 

64c 2 d 2 , 

*2r-1.2s-l - Jt 6 (2r—l) 3 (2s—l) 3 ’ r,S ~ 

Hence 


u = 


(Ac 2 d 2 


V V l,2,-l f sin 

L L ( 2 r-l) 3 ( 2 s-l ) 3 


( 2 r— \)ny . ( 2 s—l) 7 tx 
-;-sin-. 


i»= i 

Problem 2.9 Solve the equation 
d*u 1 d 2 u 
dx* + a* dt 


+ "-4 a.2 ~ 


0 < x < /, t > 0 
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where a is constant, under the conditions u = d 2 u/dx 2 = 0 , x = 0 , x = /, 
t > 0; u(x, 0) = /(x), <3u(x, t) dt = 0, t = 0,0 < x < /. 

(This can be interpreted physically as determining the transverse 
displacement of an elastic beam of length /, simply supported at its ends. 
The initial velocity of the beam is zero and its displacement is /(x).) 

Solution. The equation is not of the form of (2.2) but assuming a product 
solution X(x)1\t) gives X""-XX = 0, T" + /.a*T = 0, where A is the 
separation constant. Making the trial substitution X = Ae ,,x gives 
H* = A, i.e. n = ±A*, ±iA* hence the general solution is 

X = A cos A*x + B sin A*x + C cosh A*x + D sinh A*x. 

The conditions at x = 0 give A + C = 0 = A — C; i.e. A = C = 0; the 
conditions at x = / give 

B sin A*/ + D sinh A 1 / = 0 = D sinh /M — B sin A*/. 

Thus, if B and D are not to vanish, either sinh A*/ = 0 or sin A*/ = 0. 
The former gives A = 0 (purely imaginary values for A* are roots of 
sin A*/ = 0) which implies X = 0; hence A* = nn/l where n is any 
integer. Hence D = 0, and the appropriate representation for u is thus 

£ T(r)sin nnx/l where T"+(n 4 7 tV 7 / 4 ) T b = O.The general solution for 

n = 1 

T n is A n cosn 2 n 2 a 2 t/l 2 + B n sinn 2 n 2 a 2 t/l\ du/dt = 0 at t = 0 implies 
B n = 0 and u(0, t) = /(x) implies that A n are the Fourier coefficients of 
the sine series for f(x) in [ 0 ,f], i.e. 

2 C l .... nnx 
A „ = yj f(x)s\n — dx. 

Thus, for any given /(x), a complete series solution has been obtained. □ 

2.3 Inhomogeneous Problems It is necessary, in order to solve equation 
( 2 . 1 ) for / ^ 0 with u or du/dx taking prescribed but non-zero values on 
x = 0 , c, to modify the techniques used in the previous problems. 

In some cases / and/or the form of the boundary conditions on x = 0,c 
are sufficiently simple for a function U (often a polynomial in x) to be 
found by inspection, which satisfies both the equations and the conditions 
on x = 0,c. It is unlikely that a simple function U can be found satisfying 
whatever conditions are prescribed on y = constant also, but writing the 
unknown u as 17 + t; will give an equation of the form ( 2 . 2 ) for v with v 
satisfying homogeneous conditions on x = 0, c. The problem for v will 
be of the general type discussed earlier. 
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Homogeneous boundary conditions For problems where u or du/dx are 
prescribed to vanish for x = 0 , c (but it is not possible to use the above 
simple approach), it is only necessary to modify the previous approach 
very slightly. If the problem for/ = 0 can be solved by using a representation 
for u as a Fourier series, then a useful method to try for the inhomogeneous 
problem is to assume a Fourier series in x of the same type for f Sub¬ 
stituting both forms into the relevant partial differential equation and 
equating coefficients of trigonometric functions enables a differential 
equation to be obtained for the Fourier coefficients. Both the above 
methods are illustrated in the following problems. 

Problem 2.10 Find a function u satisfying the equation and conditions 
of Problem 2.3 except that at x = / which is replaced by du/dx = 1. 

Solution, x satisfies both the differential equation and the conditions at 
x = 0, /. Thus writing uasx + i; shows that v satisfies the same differential 
equation, and the boundary conditions on v at x = 0 and x = / are those 
imposed in Problem 2.3. The conditions u = du/dt = 0 at t = 0 give 
v = —x y dvldt = 0 at t = 0. Hence — v is the function obtained in 
Problem 2.3 and 


8 / 

M L 


n - 1 


(— iy + l ._ TtVt . .. 7TX 

(2^co S (2.-l)-s,n(2„—I)-. 


□ 


Problem 2.11 Solve d 2 u/dx 2 + d 2 u/dy 2 =2, 0<x<c, 0 < y < d, 
under the boundary conditions u — 0 on x = 0 , c; 0 < y < d\ u = 0, 
y = 0 , du/dy = 0 , y = </, 0 < x < c. 

Solution. The right hand side of the equation is sufficiently simple to 
suggest seeking elementary solutions. Particular solutions are x 2 and 
y 2 , and adding to either linear combinations of x and y will still produce 
solutions of the differential equation. If we consider x 2 + ax + h then the 
conditions at x = 0 and x = c are satisfied by b = 0, a = — c. Hence 
writing u as x 2 — cx + v shows that v is a solution of Laplace’s equation 
vanishing on x = 0 ,c; the conditions on u at y = 0 , d show that v = cx-x 2 , 
y = 0; dv/dy = 0, y = d. The appropriate representation for v is that of 
Problem 2.1, i.e. 

oc 

v = ^ [ A n cosh nny/c + B n sinh W 7 ry/c]sin nnx/c 

n = 1 

The conditions on y = d will be satisfied if A n sinh nnd/c = - B n cosh nnd/c 
and that on y = 0 shows that A n are the coefficients of the sine series of 
x(c-x) in [0,c] and thus A n have the same values as in Problem 2.1(ii). 
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Hence 


2 8c 2 V 1 

u = x -cx+— 5 - ) —- —r 

rc 3 ^ ( 2 r+l ) 3 

r = 0 




. ( 2r+l}jtx 

sin- --— 

c 

An alternative approach would be to satisfy the boundary conditions 
on y = 0 and y = d; y 2 + ay+b satisfies both the equation and these 

conditions if a = —2d and b = 0. Hence writing u as y 2 — 2 dy + w shows 

that w satisfies Laplace’s equation and the appropriate homogeneous 
conditions on y = 0, d. On x = 0 and x = c, w is equal to 2 dy—y 2 . 
The boundary value problem for w can be solved by means of a representa¬ 
tion of the form of Table l(iii). q 


Problem 2.12 Solve d 2 u/dx 2 -d 2 u/8t 2 = exp(-r), t > 0, 0 < x < l. 
under the conditions w( 0 ,f) = u(l,t) = 0 , f > 0 ; u(x,t) = du{x,t)/dt = 0 , 
t = 0,0 < x < /. 


Solution. It is not easy to find by inspection a function satisfying both 
the equation and the boundary conditions and we need to use the more 
general approach described. If the right-hand side of the equation were 

zero then the appropriate representation would be u = £ T(t) sin nnx/l 

and we therefore need to obtain a sine series representation for exp(-t) in 
[0, /]. This representation will of course be exp( - 1 ) times the sine series for 
unity, the latter series can be calculated by the methods of § 1 . 2 , and it is 
found that 

4 sin( 2 r-l) 7 rx// 

2r—1 • 0< *< l 

r = 1 


Thus substituting this series and that for u in the above differential equation 
gives 



"I . nnx 
-j2~ T n-T„ Jsin — 


_ 4 n sin(2r -1 )nx/l 
n Lt 


2r— 1 


exp(-f) 


Equating coefficients of the trigonometric functions gives 
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/2 ^2r ^2r ~ r — 0. 1> 2, • • • » 

( 2 r— 1 )V _ 4exp(-r) 

l 2 2r ~ 1 2 '~' ~ n(2r— 1 ) ’ r=1 ’ 2 ’ 


u = du/dt = 0 at t =0 imply that 7 = 7'= 0, for all n when t = 0, and 
as the equation and the conditions for T lr are both homogeneous it 
follows that T lr = 0 . The general solution for T lr _ x is 

-r _ 4/ 2 exp( — t) , , _(2r-l)7rt 

2 '~ l ~ 7r(2r-l)[(2r-l) 2 7t 2 + / 2 ] 2f -‘ / 

„ . (2r-l)7rf 

+ B 2r _ j sm - . 


The conditions at f = 0 give 
A 


4 / 2 


B 


2r ~ l 7t(2r— l)[/ 2 +(2r— l) 2 7 t 2 ]’ 
4/3 


2r ~ l 7r 2 (2r-l) 2 [/ 2 + (2r-l)V]' 
Substitution of these values into the expression for T 2r _, gives 

(2r — 1 )7rf 


4i> f i r , . 

“ " |2,-l)[P + (2r-l)V] [ e,p( - l) - 


cos- 


/ 


. ( 2 r— l) 7 tt 
r sm-— 


1 • < 
sin- 


/ 

(2r-l)7tx 

l 


□ 


n(2r — 1) 

Problem 2.13 Solve d 2 u/dx 2 — du/dt = exp( —t)sinx, 0 < x < c, t > 0, 
under the conditions u( 0 , t) = 0 , du(x, t)/dx = 0 , x = c; t > 0 , 
u(x, 0 ) = 0,0 < x < c. 


Solution. If the right hand side of the above equation were zero it follows 

from Table l(iii) that an appropriate representation for u would be 

00 

£ 7,(0 sin ( 2 n - 1 )nx/2c, and we therefore require a similar representation 

n = 1 

for sin x. This is obtained as in § 1.3, and 


sinx = 8 c cos c 


i 


(-ir 


+ i 


:sm 


(2/7— 1)7TX 


nA [(2n-l) 2 n 2 -4c 2 y~ 2c 

[It has been assumed that c is not equal to ( 2 / 7 - 1)tt/ 2 for any integer n ; 
in this case the representation of sinx would be itself.] 

Substituting both forms in the equation gives 


Z [-* 3 ^- 5 ]-^?= - *-» 

i=l L J 

® . (2n—1 )kx 

2j (2n — l) 2 Jt 2 —4c 2Sm 2c 

«i= i ' 
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and hence 


dT" {2n — 1) 2 tt 2 _ 8 ccosc(-l)"exp(-r) 

dt + 4c 2 n (2n—l) 2 7 r 2 —4c 2 

The general solution of this equation is 


2 32c 3 cos c( - 1 ) rt exp( - 1) 


T. - A.a p[-<2„— 1)>„W] + [(2 „_ |)2 ^-_ 4cST r 

and to satisfy the condition on t = 0 we require that 


32c 3 cos c(- 1)" 
[(2n — l) 2 7 t 2 —4c 2 ] 2 ' 


Hence a complete series solution for u may be written down. □ 

Problem 2.14 Solve d 2 u/dx 2 = du/dt, 0 < x < /, t > 0, under the con¬ 
ditions du(x, t)/dx = q, where q is constant, for x = 0, u(/, t) = T 0 where 
T 0 is another constant, and u(x, 0) = 0. 

Solution. Any linear function of x satisfies the equation and the particular 
choice q(x - l)+T 0 also satisfies the boundary conditions. Writing u as 
q(x — l)+T 0 + v shows that v satisfies the same equation as u, and also 
satisfies homogeneous condition for x = 0, and x = +/ and 
c(x,0) = — T 0 — q(x — l ). The appropriate representation to choose for v 
will be that of Problem 2.4 and the results can be obtained from those of 
that problem by replacing sin nx/l by -T 0 -q(x - /); this means that 


V = Yj ^n ex p[~(2« — l) 2 7t 2 t/4/ 2 ]cOS - — — 

n = 1 

where 



+ q(x - 0 ] cos( 2 n — 1 )nx/2l dx 


M-lTJo Slq_ 

( 2 n — l)n ( 2 n-l) 2 * 2 ^ 


Inhomogeneous boundary conditions When the prescribed values of 
u or du/dx on x = 0 , c are not zero and the problem cannot be reduced 
to homogeneous form, then a general, though standard, approach has 
to be used. If the corresponding homogeneous problem could be solved 
by using a particular Fourier series representation then the technique is to 
multiply the differential equation by the general member of the Fourier 
series [e.g. sin wrx/c for the case of Table l(i)] and integrate the equation 
with respect to x from x = 0 to x = c. Integration by parts of the term 
involving d 2 u/dx 2 y [for the case when Fourier series are applicable there 
will be no term involving du/dx], will lead to a differential equation for the 
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Fourier coefficients of u. This technique is illustrated in the following 
problems. 

Problem 2.15 Solve the equation of the previous problem under the 
conditions 

u( 0 ,r) = r, u(c, t) = 0 , t > 0 , u(x, 0 ) = 0 , 0 < x < c. 

Solution. The corresponding homogeneous problem [i.e. 
u{ 0,0 = u(c, 0 = 0 ] could be solved by a sine series representation and 
we thus multiply both sides of the differential equation by sin nnx/c and 
integrate with respect to x from 0 to c, i.e. 

f c d 2 u . nnx C c du . nnx 
~^r sin - dx = —sin- dx. 

Jo dx c Jo c 

Integration by parts twice of the left-hand side of this equation gives 
r du . wrx"| x = f njif n7tx"| x = f n 2 n 2 C c . nnx 

Id* c l = o c l c Jx = o c Jo c 

f c du . M 7 TX , d C c nnx 

= I — sm - dx = — u sin - dx. 

Jo dt C dt J 0 c 

00 

u(x, t) can be represented in 0 < x < c as £ &„(t)sin nnx/c, where 

# 1=1 


. 2 I*' . . . nnx , 

b = - u(x, r)sin - dx. 

C Jo c 


Thus, on applying the conditions at x = 0,c to the equation obtained 
by integration of parts. 


db_ n 2 n 2 . 2nnt 

dt c 2 c 2 

i.e. the required equation for b n . The general solution is 
4, - ^.exp(-nV«/c‘) + 2(i-^ 
and in order to satisfy u(x,0) = 0, A n = 2c 2 /n 3 n 3 giving 

u - 1 S? [ ' xp, ' n * +- l -„ sm ~■ 

n= 1 »=1 

The last series is very slowly convergent and for computational pur¬ 
poses is not satisfactory, but its exact form is given in equation (1) Appen¬ 
dix 1 and use of this result and (7) of the same appendix considerably eases 
numerical computation. In general the results of Appendix 1 are often 
useful in avoiding numerical calculation of slowly convergent series. 
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It would not have been correct in this case to have differentiated the 
sine series twice to calculate the Fourier series of 8 2 u/dx 2 as the sine series 
only represents u in the open interval and cannot be differentiated term 
by term [c.f. remarks preceding Theorem 6 ]. An alternative method would 
have been to use Theorem 6 (ii) to calculate the cosine series of du/dx and 
then applied Theorem 6 (ii) to this series to obtain the Fourier series for 
d 2 u/dx 2 . The above technique involving integration by parts is more useful 
in that it avoids the necessity for remembering formulae. 

It would also have been possible, with some ingenuity, to have reduced 
the problem for u to a homogeneous one. A simple function satisfying the 
boundary condition is u 0 = r(l — x/c), but 

du 0 = x 
8x 2 dt c 

and it is therefore necessary to add to u 0 a function /(x) such that 
f" = — 1 +x/c and with /( 0 ) = 0 = /(c); such a function is 
(x 3 — 3cx 2 + 2c 2 x)/6c. Hence writing u as 

t ^ 1 — + ^r;(x 3 — 3ex 2 + 2 c 2 x )+v 

yields a homogeneous problem for v. The form in which u is obtained in 
this case is precisely that obtained by using the results of Appendix 1 to 
transform the series solution obtained previously. □ 

Problem 2.16 Solve c 2 u/8x 2 = d 2 u/8t 2 , t > 0, 0 < x < c, subject to 
the conditions u( 0 , t) = 0 , du(x, t)/dx = sin cot, x = c, t > 0 ; 
u(x, 0 ) = du(x, t)/dt = 0 , t = 0 , 0 < x < c [to is not an odd integer 
multiple of 7t/2c], 


Solution. The appropriate representation for the homogeneous problem 
would be that of Table l(iii) and we thus multiply both sides of the equation 
by sin (2n — l)nx/2c and integrate with respect to x from x = 0 to x = c, 
i.e. 


f c 8 2 u . (2n—\)nx f c d 2 u . 

—jSin--- dx = — ^-si 

Jo 8x2 2c Jo a* 


( 2 n— \)nx . 
sm - r - dx . 


2 c 


Integration by parts twice gives 
r 8u . (2n-\)nxl x = c (2n — l)n[ 

U sm 2c 1_ 0 ~^r~[ u 


2c 

(2n— l) 2 " 2 


( 2 n-l) 7 rx] x = c 

I = o 


4c 2 


n 2 C e . (2n-l)nx J d 2 f c 

Jo 2c fc 2 J 0 


( 2 n— l)?tx , 
sm--- dx. 


2c 
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If b n is defined by u = £ 


n = 1 

prescribed conditions, that 


b n (t)sin(2n— l)nx/2c it follows, on using the 


d 2 b n (2n-l) 2 n 2 b n 2 (—l)" +1 sin cot 

~di r+ 4? = c~ 


The solution of this satisfying the conditions that dbjdt = b n = 0 at 
t = 0 is 


8 c(— 1 )" H 


(2n- \ fit 2 — 4c 2 to 2 


f . 2 toe . (2n - l) 7 tt"| 

L sma ”-p^iS sm ^r-J- 


Thus the series solution for u is obtained. 


□ 


Problem 2.17 Solve the equation 


d*u 8 4 u 8 4 u _ 

dx 4 + 8x 2 8y 2 + 8y 4 


0<x<c,0<y<d, 


under the conditions n( 0 , y) = u(c, y) = 0,(8 2 u/8x 2 ) x = o = (8 2 u/8x 2 ) xae = 
2 , 0 < y < d; u(x, 0 ) = u(x,d) = 0 , 8u/8y = 0 , y = 0 , d, 0 < x < c. 

(This problem can be interpreted physically as calculating the deflection 
of an elastic plate, two of whose edges are clamped and the other two are 
simply supported. It is also acted upon by bending moments uniformly 
distributed along the free edges.) 

Solution. This boundary value problem is of a type not previously 
considered and it can be verified that the method of separation of variables 
is not suitable. The problem can be somewhat simplified by making all 
the boundary conditions on x = 0 and x = c homogeneous and this can 
be achieved by setting u = x(x — c)+v and v will satisfy the same equation 
as u. The conditions on v show that it can be represented in [0, c] by a 

OO 

series of the form ]T V^GOsin nnx/c and, as d 2 v/dx 2 vanishes at x = 0 , c, 

n = 1 

it can be shown from Theorem 6 that it can be represented as 

oc 

— £ (n 2 n 2 /c 2 )Y n sin nnx/c for x in [0, c]. Further application of Theorem 

n = 1 

cc 

6 shows that the series representing d A v/dx 4 is £ {n*n 4 /c*)Y n sir\ nnx/c. 

n = 1 

Substituting the series in the equation and equating coefficients of 
sin nnx/c gives 

n 4 n 4 Y n -2n 2 n 2 c 2 Y;' + c 4 Y n "" = 0. 

the conditions on Y n are that Y x = 0 at y = 0 ,d and V/0) and Y n (d) are 
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the coefficients of the Fourier sine series of x(c—x). These have been ob¬ 
tained in Problem 2.1 and thus F n (0) = Y n (d) = 0 for n even and 

T 2r _,(0) = y 2r _,(d) = ^^3 , r = 1.2. 

The general solution for Y n is a linear combination of cosh nny/c, 
sinhwty/c, y cos nny/c, ysinray/c and there will be four conditions to 
determine the four constants. The calculations can be simplified by notic¬ 
ing that y„ satisfies the same condition at both boundaries, i.e. it is sym¬ 
metric about y = \d and we choose the linear combination with this 
property, i.e. 


y„ = A n cosh™(y-^d) + B n (y-\d)sinh™(y-^d). 

For n even it follows from the conditions y„(0) = Y n (d) = ^'(0) = Y„(d) = 0 
that A n = B n = 0, thus 

. (2r-l)nd d . (2r—\)nd 8c 2 

<^ Alr S j„ h ( 2Lz!M +B iSinh «tzi^ 

c 2r 1 2c 2r_1 2c 


(2r-l)nd (2r — 1 )nd n 

+ --— cosh---— = 0. 


These equations can now be solved to give A 2r _ v B 2r _ t and the series 
solution for v is then found. □ 


EXERCISES 

1. Solve 8 2 u/8x 2 + 8 2 u/8y 2 = 0, 0 < x < a, 0 < y < b, under the con¬ 
ditions u(0, y) = u(a, y) = 0, 0 < y < b; u(x, 0) = 0, du/dy = q, y = h, 
0 < x < a. where q is a constant. 

2. Obtain a solution u(p,0) of Laplace’s equation finite everywhere 
within the circle of radius a , (p,0) being the usual polar coordinates, for 
the two cases 

(i) u(a, 0) = sin 0, (ii) u(a, 9) = sin 0 ,0 ^ 6 C n 

= 0, n < 0 < 27t. 

3. Solve the equation of exercise 1 under the conditions u(0,y) = 0, 
du/dx = 0, x = a, 0 < y < b\ u(x,0) = 0, u(x,b) = 1, 0 < x < a. 

4. Solve the boundary value problem of exercise 1 when u is prescribed 
to be unity on x = 0, x = a and y = 0. 
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5. Solve d 2 u/8x 2 = du/dt, 0 < x < /, r > 0, subject to the conditions 
u(0,r) = a, u(l,t) = b, t > 0; u(x,0) = 0, 0 < x < /; a and b being 
arbitrary constants. 

6. Solve 8 2 u/dt 2 = 8 2 u/8x 2 , 0 < x < /, t > 0, subject to the conditions 
u(0, t) = 0, du/dx = 1, x = /, r > 0; u(x,0) = 0, 8u/8t = 0, r = 0, 
0 < x < l. 

7. Solve the previous exercise when the condition on x = / is u = sin t 
[nrr # / for any integer «]. 
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Chapter 3 

Fourier Integrals 


3.1 The General Fourier Integral If a function J[x) satisfies the conditions 
of Theorem 1 for all intervals then (by setting a = — c in the expressions 
of 1.1) a Fourier series representation for f(x) can be obtained in [-c,c] 
for all values of c. This series will not, however large the value chosen for 
c, represent f(x) for all x unless it has period 2c. We therefore examine 
the possibility of extending the ideas of § 1.1 to represent a function for 
all x in terms of trigonometric functions. 

In the following analysis it is convenient to use the exponential form 
of the Fourier series described immediately prior to the statement of 
Theorem 1. Any function f(x) satisfying the conditions of Theorem 1 is 
such that, except at a discontinuity, 

oo 

fix) = X 0 ,expO‘wrx/c), |x| ^ c, 



« 

3 

II 


h 


2n 

Hence 

** - h 

where 

F M iP) = f_ 


/(x)exp( — innx/c) dx 
/(x)exp( — inhx) dx. 


where ch = n. 


The ordinary definition of an integral as a limit of a sum now suggests 
that, on taking the limit h -» 0 , 




exp(i/?x)F(/}) dp 


(3.1) 


where F(p) = lim F ic \p) = H /(x)exp( - ipx) dx. (3.2) 

c -• co J - 00 

F(P) is known as the Fourier transform of/. 

The above analysis is clearly not rigorous and assumes existence of 
various limits and integrals; equations (3.1) and (3.2) are, however, valid 
under fairly general conditions, and their validity can be formalized in 
the following theorem. 


Theorem 8 If f(x ) is piecewise smooth on every finite interval of the 
x-axis and |/(x)| dx exists, then the right hand side of equation (3.1) 
represents 
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(i) /(x) at all points at which it is continuous, 

(ii) 2 [/(* r + °)+/(x r — 0 )] at any point of discontinuity x = x r . 

An important result associated with Fourier integrals is the Riemann- 
Lebesgue lemma which states that, if / satisfies the conditions of Theorem 
8 , lim F(P) = 0. 

P -* CO 

The Fourier integral is of such general applicability that tables of 
functions and their corresponding transforms have been compiled; 
some relevant results are tabulated in Appendix 2. Thus only very few 
examples of direct calculation of Fourier transforms will be given. The 
process of determining /(x) given F(P) is known as the inversion of a 
transform and this process is sometimes helped by the following theorem. 

Theorem 9 (Convolution Theorem) If / and g are functions satisfying the 
conditions of Theorem 8 and F(P) and G(P) denote their Fourier trans¬ 
forms then the inverse of FG is /(t) 0 (x-r) dt. 

Problem 3.1 Find the Fourier transform of the function f(x): J[x) = 0, 
x < 0; f(x) = exp( —x), x > 0; /(0) = Hence obtain an integral 
representation of fix) for all x. 


Solution. From equation (3.2) it follows that, since f... 0 for x < 0, 

f(W= 

/(x) satisfies the conditions of Theorem 8 , and hence 
fix) = exp iiPx)FiP) dp = J j^jj^expdPx) dp 

1 f f- U-ip)expiipx)dp . f° (1 — i/?)exp(i/?x) Jn ~| 

= 27t |_J o ^ + J 1 +/P dP \ 

P can be replaced by -/? in the second integral thus: 


Problem 3.2 Find the Fourier transform of exp(—u|x|) (a > 0) and 
hence obtain an integral representation for exp(-a|x|). 


Solution. 


FiP) = j exp(—a|x| — ipx) dx 

= (J exp( -ax-ipx) dx+ j° exp(ax-i/Sx) dx 
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D 

















1 


1 

a+ip + 


a—ip 


Thus, from Theorem 8, 


2 a 

a 2 +p 2 ' 


exp(-a|x 



exp (ipx) dp 

a 2 Tp 2 ' 


□ 


The above identity is not very easy to obtain by elementary means, but 
one method would be by means of contour integration and use of Cauchy’s 
theorem and readers familiar with complex variable theory might be 
interested in this alternative approach, [exp(ixz)/a 2 + z 2 ] is integrated 
over the closed contour in the z plane consisting of the real axis from 
z = —R to z = R and the semicircle in the upper half plane joining 
these points. There is a simple pole at z = ai with residue — i expf — ax)/2a, 
and Jordan’s lemma shows that the integral over the semicircle tends to 
zero as R -* oo, thus applying Cauchy’s theorem and taking the limit as 
R -* oo gives the result for x > 0. The corresponding calculation for 
x < 0 requires the contour to be in the lower half plane; the case x = 0 
is elementary. 


Problem 3.3 If /(x) and /'(x) are continuous functions of x satisfying 
all the conditions of Theorem 8, show that the Fourier transform of /' 
is ip F(p), where F is the Fourier transform of /. 

Solution. The required Fourier transform is 

^ f'(x)ex p( - iPx) dx = [/(x)exp( - i/Jx)] * ! * 

+ ip j ^ /(x)exp( — i/Jx) dx, [on integration by parts]. 

Since the integral of |/(x)| from x = —oo tox = oo converges it is 
necessary that /-» 0 as x -* + oo; thus the integrated terms vanish and 
the required result is proved. □ 

Problem 3.4 Show that an even function /(x) satisfying the condition 
of Theorem 8 has a Fourier representation of the form J* G(P) cos Px dp. 

Solution. The Fourier transform of /(x) is 

F(P) = /(x)exp( —i/?x) dx = 

= f; /(*)exp( - iPx) dx + f° /(x)exp( - i/3x) dx. 

x can be replaced in the second integral by -x and the evenness property 
then gives F(p) = /(x)cos Px dx. Cos px is an even function of P 
hence the same is true of F; a particular example of this general result has 
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been derived in Problem 3.2. /(x) is such that Theorem 8 can be applied 
and thus 


m 


i_ r® 

~ 2?r 

J — * 

_L f* 
= ^Jo 


F(p)exp(ipx) dp 


F(/?)exp(i/?x) dP+ 


1 f° 
2zr J _, 


F(P)exp(iPx) dp. 


Replacing p by — ft in the second integral and using the evenness property 
of F(p) gives f(x) = (\/n)$™ F(fI)cos Px df}, which is of the required 
form. □ 

It can similarly be shown that, for / an odd function of x, 

/(x) = (1/7r) j 10 G{P)s'm Px dp where G{P) = 2 l°° /(x)sin^xdx. 

Jo Jo 


3.2 Fourier Cosine and Sine Transforms By constructing even and odd 
extensions of a function /(x) defined for x > 0, one can (from Problem 
3.4) obtain integral representations involving only sines or cosines [this 
is analogous to obtaining sine and cosine series representations on a 
finite interval]. Hence we have the following theorem. 

Theorem 10 For a function /(x) satisfying the conditions of Theorem 8 
for x > 0 the Fourier cosine and sine transforms F c (P) and F s (P) defined by 


F c (fi) = } ” /(x)c° s Px dx, (3.3) 

F s (p) = j* /(x)sin Px dx (3.4) 

are such that at all points x > 0 at which / is continuous, 

/(x) = (2 /tt) j® F c (0 )cos px dp, (3.5) 

f(x) = (2/rr) j* F s (P)sinPxdp. (3.6) 


At those points x = x r at which / is discontinuous the left hand sides of 
equation (3.5) and (3.6) have to be replaced by i[/(x r +0)+/(x r -0)]. At 
x = 0 the right-hand side of equation (3.5) converges to /(0 + 0) and that 
of equation (3.6) to zero. 

Problem 3.5 Calculate the Fourier cosine and sine transforms of 
exp(-x)cosx. 

Solution. 

F c (P) = J® exp( —x)cosxcos px dx, F s {P) = J* exp( - x)cos x sin Px dx. 
The most direct method of calculating these integrals is by expressing 
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cos fix and sin fix in terms of exp( ± ipx) and taking cos x as the real part 
of expfix). Thus, on defining G(P) by 

1 1 + i(P+l) 


G(P) 


-i; 


exp( — x + /x + ifix) dx = 


we have that 

Ftf) = *Re[G(/?) + G(-/?)], 
Hence 


i-np+i) 2+ip+p 2 ’ 
F s (p) = iRe -i[G(/?)-G (-PM 
P 2 + 2 


1 Re 2( 1 4- Q(2 4- p 2 ) —4f/? 2 
' {P 2 (2 + p 2 ) 2 — 4p 2 p* + 4' 

lRc4i(\ + i)p + 2P(2 + p 2 ) P 3 

sW 2 p*+4 P*+4 

From these expressions and Theorem 10 are obtained the identities 

*• (p 2 +2) 


exp( —x)cosx 
exp(-x)sinx 


- 2 -[ 

_ 2 f® p 3 si 
■-Jo 


o ^+4 
sin px 
+4~ 


cos Px dp 
dp. 


f(P)cosPxdp = 


Problem 3.6 Solve the integral equation 

1, 0 < x < 1, 
0, x > 1. 

Solution. Theorem 10 effectively states that if 
[”/(/J)cos Pxdp = F(x) 


£ 


then J (P) = (2/tc) j* F(x)cos Px dx, 

and applying this result gives 

2 f 1 2sin/? 

JW = ~ cos px dx = — -2-, 

- Jo *0 

Substituting this into the original equation and setting x = 0 shows that 

f® sin P . 

i.- 0 -** ° 

Problem 3.7 Find the Fourier cosine transform of exp(—x 2 ). 

Solution. We have F c (P) = J® cos px exp( - x 2 ) dx. The simplest method 
of evaluating this integral is by a complex variable method but we shall 
consider an alternative method which uses more elementary ideas. The 
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integral is uniformly convergent for x > 0 and thus dFJdp can be evalu¬ 
ated by differentiating under the integral sign, giving 

dFJdp = - j.* x sin Px exp( —x 2 ) dx = \ sin pxd exp(-x 2 ) 

= -Ws 

(on integration by parts). 

Hence F = A exp (-j/? 2 ) where A is a constant. A may be found by 
setting P = 0 and using the known result 2 J* exp(-x 2 ) dx = Jit. The 
final result is 2F C (P) = Jit exp (—jp 2 ). □ 

[For readers familiar with complex variable theory the following 
direct method might be of interest. As e~ xl is an even function of x it 
follows, by an analysis exactly similar to that of Problem 3.4, that 

2 F(P) = P° exp( — ipx — x 2 ) dx = exp(- tf 2 ) f°° exp( - x 4-|i/?) 2 dx. 

c J - CO J-00 

exp(—r 2 ) is an analytic function of z and thus, integrating it round the 
rectangular region with vertices at z = ±R, ±R+^iP, gives 

exp (~z 2 )dz+ + *'* ex p( _ z 2 )dz + * exp(-z 2 )dz-t- 

4- f R exp( — z 2 ) dz = 0. 

J-K+4V 

As R -* oo the second and fourth integrals tend to zero, and setting 
z = x+jiP in the third integral gives 

p exp (~z 2 )dz= j* exp[-(x4-§//?) 2 ] dx. 

The left-hand side of this identity is known to be Jit and thus F c {p) can 
be calculated.] Q 

Problem 3.8 If /, f and /" are continuous functions of x satisfying all 
the conditions of Theorem 10 determine the cosine and sine transforms 
of f" in terms of the corresponding transforms of / and the values of / 
and f at x = 0. 

Solution. The cosine transform of /" is G c (P) defined by 
GAP) = j" /"cos Px dx 
and integrating by parts we get 

G C (P) = [f cos px] x x lo+P j* /'sin px dx 

= [/'cos PxY x i;+PU sin PxY x :^-P 2 \ * / cos Px dx 
(on a further integration by parts). /, /' both vanish as x -* oc and thus 
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G c (/J) = -[/'] x = 0 -/f 2 F c (/J). Repeated integration by parts of the cor¬ 
responding integral for G S (P), the sine transform of /, gives 

gw = Pin xm0 -P 2 F S (P). □ 

3.3 Solution of Boundary-Value Problems Using Fourier Transforms 

Fourier cosine and sine transforms are useful for the solution of boundary 
value problems for 

d 2 u d 2 u du 

a i-j^ +a 3 u+b i(y)fy2 +b 2(y)fy+ b 3iy)u = /(*,y), x > 0,0 < y < d, 

where a v a 3 are constant, with u-*0asx-*x and either u or du/dx 
prescribed on x = 0. The general method is to take either the sine or 
cosine transform of the above equation; this gives (using the results of 
Problem 3.8) 

a l [-p 2 U c -(du/dx) x . 0 ]+a 3 U c + j *Lu cos Px dx = F e , (3.7) 

a i[~P 2U s + P( u ) x = o] + a 3 U s+ \*LusinPxdx = F s , (3.8) 

where capitals denote Fourier transforms, the suffixes c and s have the 
same meaning as before, and L denotes the operator b t d 2 /dy 2 + b,d/dy+b 3 . 
In deriving the above results it has been assumed that u, du/dx, d 2 u/dx 2 are 
continuous in x and satisfy the conditions of Theorem 8. It is also assumed 
that differentiation with respect to y may be interchanged with the integra¬ 
tion operation in equations (3.7) and (3.8) giving 

a xl-P 2 V~ 0K0XUJ+O, U c +LU e = F c , (3.9) 

^l~P 2 G s +P(u) xm0 ]+a 3 U,+LU S = F s . (3.10) 

Thus, if du/dx is known at x = 0, equation (3.9) becomes a differential 
equation for U c regarded as a function of y; similarly, if u is prescribed 
at x = 0, equation (3.10) gives the differential equation satisfied by U s . 
Once either U c or U s has been found then Theorem 10 can be used to 
find u. Various assumptions have had to be made about the unknown 
function and it has to be verified that the function eventually determined 
satisfies these conditions. The above general procedure is formally 
identical with that described in §2.3 for the solution of inhomogeneous 
boundary value problems using Fourier series. 

If x ranges from — x to x and the conditions on u are u -> 0 as x ± oo, 
then the exponential Fourier transform has to be used. In this case it is 
possible to consider the more general equation 

d 2 u du 

a ifo2 + a 2fa + a 3 U + Lu = f, (3.11) 


48 



where a v a 2 , a 3 are constant. Taking the Fourier transform of equation 
(3.11) gives, using the result of Problem 3.3, 

— p 2 a 2 U + iPa 2 U + a 3 U+LU = F. (3.12) 

Various assumptions concerning continuity and interchanging operators 
have again been made in deriving equation (3.12). This is a differential 
equation for U and, after it has been solved, u may be found from Theorem 
8 . 

Problem 3.9 Solve d 2 u/dx 2 = du/dt, -oo<x<oc, t>0 under the 
conditions u -» 0 as |x| -* oo, u(x,0) = exp(-x 2 ). 

Solution. The appropriate transform to use is clearly the exponential 
one, and taking the transform of the equation gives —p 2 (U) = 
dU/dt. At t = 0, u = exp(-x 2 ) and hence U = exp(-i0x-x 2 )dx 
at t = 0, and it follows from Appendix 2, equation (5) that U = 7t*exp( -\P 2 ) 
when t = 0. Solving the differential equation for U gives 

U = A(P)txp(—P 2 t) 

where A can be a function of p. The condition on r = 0 shows that A(p) 
is equal to the value of U at r = 0 and hence 

U = ;t*exp[-^/? 2 (l + 4t)]. 

The inverse of 7r i exp(— \P 2 ) is exp(—x 2 ) [Appendix 2, (5)] and from Appen¬ 
dix 2, (1) it follows that the inverse of exp(—|A 2 /? 2 ) is x -1 exp( —A“ 2 x 2 ). 
Hence setting A = (l+4t)* shows that u = (l+4r)"*exp(-x 2 )/(l-l-4r). 
It can be verified in this case, by direct substitution, that u is the required 
solution. □ 

Problem 3.10 Solve d 2 u/dx 2 +d 2 u/dy 2 =0, 0<x<oo, 0 < y < d, 
under the conditions u(x,0) = u(x,d) = 0, x > 0; u(0,y) = l—y/d, 
0<y<d;u-*0asx-*oo. 

Solution. This problem could be solved by either representing u as a 
Fourier sine series in y or taking the Fourier sine transform with respect to 
x, and we choose the latter approach. Taking the Fourier sine transform 
of the above equation gives 

U s = -P(l-y/d), [c.f. equation (3.10)]. 

The general solution of this equation is 

U s = A cosh Py + B sinh Pv+(d-y)/pd. 
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and the conditions U s = 0 for y = 0 and y = d [which follow from 
u = 0 on y = 0, rf] yield 


Hence 


"Jo 


v _ Sinh Piy-d) (d-y) 
5 ft sinh fid pd 
sinh P(y - t/)sin Pxdp 2 (d-y) 
P sinh pd + n d 


r 


sin Pxdp 


jo P ' 

Setting a = y-d,y = d in Appendix 2, equation (8) gives 

f* sinh x(y-</)sin px ^ f. [n, „~| . fnp\\ 

J. — sinhJ* dx - a " | an L25 0 ' _ ‘ 0 J* anh l25;j 

and interchanging the variables p and x gives the first term in u. The 
second term in u can be evaluated by using the identity derived at the end 
of Problem 3.6 (replacing p by Px) and thus finally 


i y 2 

u = 1— - — tan 


■H~s} 


Clearly this function satisfies all the conditions imposed. 


Problem 3.11 Solve the equation of the previous problem in 0 < y < d, 
x > 0, subject to the conditions u(x,0) = 0, u(x,d) =1; x > 0; 
m( 0, y) = 0, 0 < y < d. 


Solution. The sine transform is still applicable and now 


d 2 U , 
6y 2 


-p 2 u s 


= 0. 


As u vanishes on y = 0, U s must also vanish on y = 0; u(x,d) = 1 and 
hence U s on y = d is the sine transform of unity. However, unity does not 
satisfy the conditions of Theorem 10 and it does not seem possible to 
calculate its Fourier sine transform. It follows however from the last 
identity of Problem 3.6 by replacing p by Px, (x > 0), that there exists an 
F(P) (in this case 2 /np) such that J® [F(0)sin px]dp = 1 . F(P) is not abso¬ 
lutely integrable and thus Theorem 10 cannot be used to invert the above 
result. We therefore seek an integral representation for u of the form 
u = Jo G(P,y) sin Px dp. Substituting this in the equation for u shows that 
G satisfies the same equation as U s above. The boundary conditions give 
G(P,d) = 2/nP and G(p.o) = 0, and hence G = 2 sinh Py/p sinh pd and 

u _ 2 f* sinh Py sin px dp 
n J o P sinh pd 

This integral is of the same form as that of Problem 3.10 with y-d now 
replaced by y and hence 
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□ 


2 -i 

u = - tan * 

7T 

(For all problems with homogeneous boundary conditions the above 
method of assuming an appropriate integral representation of the solution 
can be used instead of taking the transform of the equation. The two 
methods are entirely equivalent but, provided the relevant transforms 
exist, taking the transform of the equation has the advantage of also being 
applicable for non-homogeneous boundary conditions. This is also the 
method used for solving ordinary differential equations using Laplace 
transforms.) 

Problem 3.12 Solve the equation of Problem 3.9 for x > 0, t > 0 
subject to the conditions « -* 0 as x -» oo, n(x,0) = 0 and du/dx = 1, 
x = 0, r > 0. 

Solution. In this case the only possibility is the cosine transform and 
—P 2 U c — 1 = dUJdt [c.f. equation (3.9)]. The general solution is 

U c = -^ + A(p)cxp-P 2 t 

and, as U c = 0 for t = 0, it is necessary that A — /? -2 . Thus applying 
Theorem 10 gives 

2 f* cos/?x(exp(-/? 2 r)— 1) 

“";Jo - w - * 

Appendix 2, equation (12) gives (interchanging x and /?) 

f™ exp(-a/? 2 )cos pxdp = -(7r/a) i exp( — x 2 /4a) 

and integrating both sides of this with respect to a from a = 0 to a = t 
gives 

p-y’W f'«-*exp(-x>/4a)fc 
Jo P Jo 

Making the change of variable v 2 = x 2 /4a enables the integral on the 
right-hand side to be rewritten as 

-x^ x 2t Qxp(-v 2 )/v 2 dv = 2r*exp( —x 2 /4r) + 2x [* exp( — v 2 )di\ 
on integration by parts. 

Hence u — — 2r*exp(— x 2 /4t)/n* + x Erfc(x/2f*), 

where Erfcz, the complementary error function, is defined by 
Erfcz = 27r“* J* exp(-r 2 )dr. 




□ 
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Problem 3.13 Solve 

3m £ 2 m 

0<r< a, z>0 
under the conditions w(0,z) finite, u(a,z) = 0 and u(r y 0) = 1. 


5 2 u 1 8u d 2 u _ 

dr 2 + r dr + dz 2 


Solution. This problem can be solved by taking the sine transform with 
respect to z and this gives 


8 2 U S 

dr 2 




A particular integral of the equation is 1/j? and the corresponding homo¬ 
geneous equation is Bessel’s equation of order zero in the variable i/?r. 
The appropriate finite solution of this is J 0 (if}r) which is normally written 
as / 0 (j8r). 

The solution satisfying the condition on r = a is thus 


U 


5 



Wl 



The first term is equal to unity (c.f. last identity of Problem 3.6) and the 
second integral can be transformed by use of Cauchy’s theorem to yield 
the form of solution obtained by using a Fourier-Bessel expansion 
(cf. Chapter 5). From the numerical point of view the infinite integral 
form is as convenient as the infinite series form. □ 

Problem 3.14 Solve d 2 u/8x 2 + d 2 u/dy 2 = 0, y > 0, all x, under the con¬ 
ditions w(x,0) = /(x), u -+ 0 as x 2 +y 2 -+ oo, where / is a function 
satisfying the conditions of Theorem 8. 

Solution. The two possible methods of solution are the exponential 
Fourier transform with respect to x or a sine transform with respect to y. 
The latter will lead to an inhomogeneous equation for the Fourier trans¬ 
form and this generally presents more technical difficulties. We therefore 
take the exponential Fourier transform with respect to x and have 
d 2 U/dy 2 -(} 2 U = 0. In order to satisfy the conditions at infinity 
U = exp( — \/}\y) and to satisfy the condition on y = 0, A(f}) must be 
equal to F(/?), the Fourier transform of f. Hence U = F(/?)exp( - \p\y). 

Theorem 8 now gives a formal expression for u but a simpler one is 
obtained by applying Theorem 9 (the convolution theorem) and using 
the result that the inverse of exp( — \P\y) is y/n(x 2 -\-y 2 ) [Appendix 2. (4)]. 
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Hence 




mdt 

(x-£) 2 +y 


2 • 


Problem 3.15 Solve 


d\ _ 

dx 4+ dx 2 3y 


8*u d*u 

3-V’ 0 ’ 


y > 0, all x, 


under the conditions du/dy = 0, y = 0; (du/dx) y = 0 = /(x), where /(x) 
is a function satisfying the conditions of Theorem 8. Also u 0 as x, 
y -> oo. 

Solution. This equation is not of the general type that we have considered, 
but by analogy with our previous examples it would seem reasonable to 
attempt a Fourier transform and use of the exponential form avoids any 
complications with inhomogeneous terms. Thus 

o.rj 2B 2 S 2 U 0*17 

P U ~ 2P ~df ' 

the general form of solution tending to zero as y -* oo is 
U = [/4 + By]exp( — |/?|y), the condition on du/dy on y = 0 gives 
\P\A = B. The condition on du/dx on y = 0 gives ipU = F, where F is 
the Fourier transform of / and thus 

r 

U- + \p\y)sxp(-\p\y)F(p). 

The inverse of the term multiplying F(P) is 

0+1 P\y). 


-fj; 


p 


exp(-\P\y)exp(iPx) dp 
sin px 


-if 

n Jc 


- r si 

n Jo 


p exp (~Py)dP+ z I sin Px exp( — Py) dp. 

From Appendix 2(9) the first term of the above equation is equal to 
(l/ 7 r)tan - { (x/y) and Appendix 2(12) shows that the second integral is 
equal to x/n(x 2 + y 2 ). The solution for u can now be written down by 
means of theorem 9. □ 


Problem 3.16 Solve Problem 3.12 when the condition on x = 0 is 
du/dx = hu where h is a constant and u(x 0) = exp( - x). 

Solution. If h were zero then the appropriate transform would be a 
Fourier cosine transform and for h -* oo the Fourier sine transform 
would be appropriate but neither is suitable for other values of h. It 
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therefore seems worth considering the possibility of using a combination 
of sine and cosine transforms and we consider U* defined by 

U* = a usin fixdx + b (* u cos fix dx, 

where a and b are functions of p. To obtain the differential equation 
satisfied by U* it is necessary to calculate 

f® . d 2 u T 

J (a sin fix+b cos fix) dx = j^(a sin px + b cos Px) —J 


= j^(a si 


-j: 


(Pa cos px—pb sin Px) ^ dx 
ox 


. du 


sin Px+b cos Px) — - (pa cos Px—pb sin Px)u 


■r 


-p 


j: 


(a sin Px + b cos f}x)u dx. 

0 as x -♦ oc; 


f 


For any of the integrals to converge it is necessary that u 
this means the integrated terms vanish as x -* oc and 
x> ^2^ r g ”i 

(a sin Px+bcos Px)j^dx = -b^ + aPu -p 2 U*, 

and, in view of the boundary condition on u at x - 0, the first term 
vanishes if aP = hb and the simplest choice of a and b gives a = h; 
b = p. Multiplying the equation for u by h sin Px+P cos Px and integrating 
with respect to x from 0 to oc formally gives dU*/dt = ~p 2 U*. 

At t = 0, U is exp(-x)and hence 

U* = j* (h sin Px + P cos /Jx)exp( - x) dx = P(h+l)/(l +0 2 ), 

Thus U* = (l+h)pcxp(-p 2 t)/(l+p 2 ). It is now necessary to determine 
the inversion formula for this transform, we have 

U* — J sin Px dx + J u—(sin Px) dx = J ^hu-^jsin Px dx, 

(on integrating the second term by parts). From this expression it is seen 
that U* is the sine transform of hu — du/dx and hence 


8x * Jo 


U* sin px dp. 


This equation can be rewritten as 

— ~(ue~ hx ) = -e~ hx 
ox n 


r 


U*sinpxdp 
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and replacing x by a variable w and integrating both sides with respect 
to w from w = x to w = oc gives 


_ 2 p* (h si: 

U ~ * Jo 


sin Px + P cos px)U* 


h 2 +p 2 


dp. 


The above derivation of the inversion formula has made various assump¬ 
tions concerning interchanging orders of integration, but it can be 
verified that the inversion formula is valid for all functions /(x) satisfying 
the conditions of Theorem 8. Thus the solution of the present problem 
can be obtained in the form of an infinite integral. □ 


EXERCISES 

1. If /(x) is exp(2x) for x ^ 0 and zero for x > 0 show that 

2 cos Px - P sin Px jp 


fix) 


i f" 

~ 71 Jo 


P 2 + 4 


2. Obtain the Fourier sine transform of exp( - x)sin x, (x > 0) and hence 
show that,for integer nj” (x sin njrx)/(x 4 + 4) dx — 0. 

3. Solve the integral equation 

= 0 , 0 ^ p ^ 1 , 

f* /(x)sin Pxdx =1, 1 < P ^ 2, 

*'0 

= 0 , P >2. 

4. Solve du/dt = d 2 u/dx 2 , x, t > 0, subject to the conditions 
u(0, t) = sin cot, u -» 0 as x -» oo, t > 0; u(x,0) = 0, x > 0. 

5. Solve d 2 u/8x 2 -3u = du/dt, x, t > 0 under the conditions u(x,0) = 0, 
x > 0; du/dx = e~\ x = 0, t > 0; lim u(x,t) = 0, t > 0. 
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Chapter 4 

Generalized Series Expansions 


4.1 Introduction In Chapter 1 it was shown that a wide class of functions 
could be represented by series of trigonometric functions and we now 
examine the extension of this principle to represent a function as an 
infinite series of certain other types of functions. In practice it is only 

useful to consider a representation of / on [a, b ] by a series of the form 
00 

£ a n (/) n where the satisfy a relation of the form 
» = 1 

r^n^m dx = 0 ’ # M, (4.1) 

J a 

where p is a given function. In this case the functions </> n are said to be 
orthogonal with weight p on \a,b]. All the Fourier series previously 
considered are of this form with p = I (c.f. Problem 1.3). 

If it is assumed that, for a given set of functions <f > n , a function / is such 
that 

/ = X a n^ a < x <b, (4.2) 

n = 1 

then, multiplying both sides of this equation by p<p m and integrating with 
respect to x from x = a to x = b gives 

\ h pf<t> m dx = £ af p<t> H <t> m dx, 

J a w — 1 Jq 

(assuming that the order of summation and integration may be inter¬ 
changed). Equation (4.1) now shows that 

£ pf K dx = f pK dx > (4 - 3) 

thus yielding an expression for a m . The above analysis is purely formal as 
it need not be true that, for a given set of equation (4.2) should hold, 
for any /. (If the <j) n were chosen to be l,cos 2x,cos3x,...on [0,7i] then 
for / = cos x, equation (4.3) would give a m = 0, thus yielding an obvious 
contradiction and showing that the above set of functions would not be 
suitable to use in representing a function on [0, n].) In those cases where 
equation (4.2) is meaningful the coefficients defined by equation (4.3) are 
known as the generalised Fourier coefficients of / and the resulting series 
representation as a generalized Fourier representation of /. The formal 
procedure given in Problem 1.3 is a special case of the method outlined 
above. 
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The above example shows that a necessary condition for a representation 
of the form of the right-hand side of equation (4.3) to be valid is that there 
exists no function orthogonal to all the $ n . In making statements of this 
nature concerning the existence of a function, one has to be somewhat more 
precise and specify the class of functions considered (e.g. functions con¬ 
tinuous on an interval, or the class of integrable or square integrable 
functions over an interval). The class of functions is sometimes said to 
constitute a function space. If the functions </> n are orthogonal with weight 
p on [a, fc] and such that there is no function / of the same class, ortho¬ 
gonal to all the (j) n and such that £ pf 2 dx # 0, then the set of functions 
(j) n is said to be complete within the given class. If the series of equation 
(4.2) with a m defined by equation (4.3), converges to /, except at possibly 
a finite number of points in [a, ft], for all functions / of the class con¬ 
sidered, then the set (f> n is said to be closed in the sense of pointwise 
convergence. This terminology is introduced here to make the reader 
aware of it but the concepts will not be used in the text. 

We now examine briefly two general classes of problems which lead 
naturally to the idea of expanding a function in terms of two particular 
sets of functions. We consider first the problem of solving Laplace’s 
equation in a spherical region when axi-symmetric boundary conditions 
are prescribed on one, or two, spherical surfaces. Laplace’s equation 
becomes, on assuming axial symmetry, 


1 d f 2 1 d 

r*dr\ r d?) + 71in0de 



= 0, 


where r, 6 are the usual spherical polar coordinates. Writing u as F,(r)F 2 (0) 
and using the procedure of separation of variables gives 



and, on writing x as cos 6, equation (4.4) becomes 


d_ 

dx 



+ XF 2 


= 0. 


(4.5) 


In physical problems F 2 has to be finite for all 0 (i.e. |x| < 1) and it can be 
shown from the theory of ordinary differential equations that equation 
(4.5) only has such finite solutions when k = n(n +1), where n is an integer. 
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In this case there is a finite solution, the Legendre polynomial, P n (x), of 
degree n satisfying 


d_ 

dx 


[<->£] 


+ n(n+ 1)P = 0, 


(4.6) 


and the condition P n (l) = 1. 

When 2 = n(n +1) independent solutions for F, are r" and r ~ n ~hence 

ao 

an appropriate series representation for u is £ + B r r~ H ~ l )P H (cos 0). 

m = 0 

Thus if u is prescribed for r = a and r = b in terms of 9 and if any function 
of 9 can be expressed as a series of Legendre polynomials then A n and 
B n can be found. Boundary value problems of this type lead to the concept 
of expanding a function as a series of Legendre polynomials. The formal 
properties of such an expansion are given in a subsequent section but it is 
convenient to summarise some of the simpler properties of these functions 
at this stage. 

An explicit form for P n (x) is provided by Rodrigues' formula 


W = 


1 d” 
2 Vdx” 


(x 2 -ir, 


(4.7) 


some particular cases of this being 

P 0 (x) = 1, (4.8) 

Fj(x) = x, (4.9) 

2P 2 (x) = 3x 2 —1. (4.10) 

Also 


P„(-x) =(-l ) n P n (x), (4.11) 

P 2" (0)= * 2^.4 2 ^2n ' ’ n = 1 ’ 2 ’-’ (412) 


P 2(I+1 (0) = 0, (4.13) 

P' 2n ( 0) = 0, (4.14) 

P' n+1 (0) = (2n + l)P 2 „(0), (4.15) 

P o (0) = P',(0) =1, (4.16) 

(1 - 2xf + f 2 )~ * = J f"P„(x). (4.17) 

n = 0 


The Legendre polynomials are orthogonal with unit weight on [ — 1, 1] 
and 



2 

(2h + 1) 


(4.18) 
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Laplace’s equation in cylindrical polar coordinates (p, <j>, z) becomes, 
assuming axial symmetry, 

1 d ( du\ d 2 u n 
pdp\ dp) + dz 2 

writing u = G,(p)G 2 (z) and using separation of variables gives 




(4.19) 


Equation (4.19) can be rewritten as 

dy y dy 1 

where y = this is Bessel’s equation of order zero with solutions 
J Q (y) and Y 0 (y). For finite solutions in boundary value pFo6Iems~Tor 
regions including the origin Y 0 must be excluded. For the particular 
class of problem when u is required to vanish on p = 1 it is necessary 
that / 0 (A*) = 0 and thus A* must be one of the zeros of / 0 (y). There are 
an infinite number of these which we shall denote by j 0n . The equation for 
G 2 is soluble in simple terms and thus the appropriate series representation 

X 

is X M n cosh j 0n z+B n sinh/ 0 „ z)J 0 (j 0n p). If u is prescribed on z = a 

n = 0 

and z = b and if any function of p can be expanded as a series of J 0 (j 0n p) 
then A n and B n may be found. 

In §3.3 the general concept of expanding a function in [0,1] as a series 
of terms of the form J v {j vn x\ where J y {j vn ) = 0, is examined in detail. Some 
of the simpler properties of Bessel functions are summarised below. 

The Bessel function 7 v (Ax) satisfies the equation 

- °- < 4 - 2 °» 


and 

Also 


7 v Ux) = I 


v + 2r 


(-mix) 

o r\(v+r)\ 


(4.21) 

(4.22) 

(4.23) 

(4.24) 

= 2J M), (4.25) 

where the dash denotes the derivative with respect to the argument. The 
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d[xV,(Ax)]/dx= AxV v _ ,(Ax), 
d[x-V s (Ax)]/dx = — 2 x - V v+ ,(/jc), 

A_,(x)+y, + 1 (x) = -/.(X), 


E 




















set of functions Mj x „x) [/? = 1 , 2 ,...] are orthogonal with weight x on 
[ 0 , 1 ] and 

f l x/ v0'. n x)dx = K 2 + 1 (/„)• (4.26) 

Series involving J v (j vn x) are generally known as Fourier-Bessel series. 


4.2 Expansion in Terms of Legendre Functions The basic result is given 
by the following theorem. 

Theorem 11 If / is piecewise smooth in [ — 1,1] then the series 
Y. a „ P„(x),where 2a n = (2n + 1) Jlj /(x)P n (x) dx, converges for |x| < 1 

n = 0 

to f(x) at all points at which / is continuous and to 2 [/(x + 0 )+/(x— 0 )] 
at points of discontinuity. At the end points x = +1, x = -1 the series 
converges to /(1 — 0) and /(—1+0) respectively. The particular form 
occurring for a n can be deduced directly from equations ( 4 . 3 ) and (4.18). 

J f(x)P n (x)dx = f(x)-^- n (x 2 -l) n dx, [from equation 4.7] 

i r d n ~ i , i 1 if 1 


1M- 1 


df_d^ _ 

dxdx n ~ 1 


(x 2 — 1 )" dx , [on integration by parts.] 


(x 2 - l) n and its first (n- 1 ) derivatives vanish at x = 1 , x = - 1 and thus 
the integrated part in the above expression vanishes. The process can be 
repeated n times giving 

£ J‘ ( £(* 2 -l Tdx. (4.27) 

Equation (4.27) is often useful in calculating coefficients in Legendre 
expansions. 

Problem 4.1 Find the Legendre expansion on [ — 1,1] of x 2 . 

Solution. Equation (4.27) and Theorem 11 show that the coefficient 
of P„(x) is given by 

(-1)"i 


a _ = 


l)"( 2 n + 1 ) f 1 d n 

2 n+, n! J_, dx” 


(x 2 Xx 2 -iy dx. 


The nth derivative of x 2 is zero for n > 2 hence a n = 0. n > 2 and 
a i = I flj (x 2 - 1) 2 dx = §, a, = -f J‘ 4 x(x 2 - 1 ) dx = 0 , 


“•"tJ’, 


x 2 dx = 3 . 
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Hence x 2 = §P 2 +jP 0 ; this particular result could have been obtained 
more directly from equations (4.8) and (4.10). □ 

It can be shown similarly that the coefficients of P n in the expansion of 
x m vanish for n > m. In this case the Legendre expansion would only have 
a finite number of terms and is best obtained by comparing the coefficients 
of powers of x in the expansion x m = £ a n P n (x), where the explicit 

n = O 

expressions for the Legendre polynomials are substituted on the right- 
hand side. 


Problem 4.2 Expand cos 20 as a series in P n (cos 0) in 0 ^ 0 ^ n. 

Solution. In this example x has been replaced by cos0 and the same 
substitution should be made in the integral for a n . It is, however, simpler 
in this case to use cos 20 = 2cos 2 0-l and equations (4.8) and (4.10). 
Thus 

cos 20 = 2cos 2 0-l = f[2P 2 (cos0)+1] — 1 
= fP 2 (cos<£)-^ = fP 2 (cos0)-’P o (cos0) (since P 0 =1). □ 

Problem 4.3 Find the Legendre expansion in [— 1,1] of/(x): 

/(x) = 0 , (-1 < x < 0 ), /(x) = 1 , (0 < x ^ 1 ). 


Solution. In this case a n = $(2n + 1)Jq P n (x)dx. The integral can be 
evaluated by using equation (4.7) but a more direct approach is to use 
equation (4.6). giving 


a 


n 


(2n + 1 ) 
2n(n + 1 ) 



dx 


2n+ 1 
2n(n + 1 ) 


p;,( 0 ). 


PJ0) can be calculated from equations (4.13) to (4.16), which show that it 
is zero for n even (n # 0 ) and 


(-0^ + 3 1.3.5 — (2/c —1) 

a2k+i (4k + 2X2k + 2) 2.4. 6 ...2 k ,A ‘“ 

°l = 4> fl O = 2- 

By theorem 11 the series should converge to ~ at x = 0, which is a point of 
discontinuity and the above results confirm this as, from equation (4.13), 
P 2n+ ,( 0 ) = 0 and the value of the series at x = 0 is a 0 = □ 

Problem 4.4 Show that it is possible to represent a function /(x) de¬ 
fined in [0,1] as a series of Legendre polynomials P n (x) of even order. 

Solution. Since Theorem 11 refers to the interval [ — 1,1] it is necessary 
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to find a function F(x) equal to fix) in [0,1] and such that the coefficients 
of odd order in the Legendre expansion of F are zero, i.e. 

\' F(x)P n (x)dx = j ‘ P„(x)F(x)dx + j° ] P n (x)F(x) dx 

= £ P n (x)\_F{x )+(-1 )"F( - x)] xdx = 0, n odd 


on replacing x by -x in the second integral and using equation (4.11). 
Thus F( — x) = F(x) and hence F(x) has to be defined for - 1 < x < 0 
as/(—x). The coefficients a ln of P 2n (x) are given by 

4n + l f 1 f 1 

a 2n = —J~ J F (*) P 2.( X ) dx = 4/1+1 I /(x)P 2b (x) dx. 

The series converges at x = 0 to ^[F(0 + 0)+F(0 - 0)] = /(0). □ 

Similarly fix) can be represented as a series of Legendre polynomials 
of odd order by defining F(x) = -f(-x) for -1 =$ x < 0. 

In this case the resulting series converges at x = 0 to 
i[F(0+0) + F(0-0)] = 0. 

Problem 4.5 Find the series representing x on [0,1] in terms of Legendre 
polynomials of even order. 

Solution. This is a special case of the previous example with f(x) = x 
and hence a 2n = (4« + 1) fg x P 2n (x) dx. The integral is best calculated by 
using equation (4.7) giving 


a 2„ = 


4/t + 

2 2 "(2/i) 


1 f 1 d 2 " , , 

^Jo X ^ (X ' 


l) 2 " dx 


(4/i+D rr </ 2 "-'(x 2 —d 2 "T =i f 1 d 2n ~ i (x 2 — i) 2 " ) 

2 2 "(2/i)! ||_ X dx 2 - 1 1. 0 J 0 /fx 2 " -1 dx \ 

n # 0. 


The integrated part vanishes at both limits and thus, after a further inte¬ 
gration. 


(4/i+l) 

2 2 "(2/i)! 


d 2 "" 2 

dx 2 "- 2 



n ± 0, 


the expression in square brackets vanishes at x = 1 and can be evaluated 
at x = 0 by expanding (x 2 - l) 2 " by means of the binomial theorem, yield¬ 
ing 


4/i+l d 2 " -2 ^ x 2r (— l) 2 " -r (2/i)! 
' 2 " “ 2 2 b (2/i)! dx 2 "' 2 o r!(2/i —r)! 


x = 0. 
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It is only the term x 2 " 2 that produces a non-zero contribution after 
differentiating 2zi—2 times and setting x = 0, and thus 

(4/i+ 1X2/1 —2)!(— 1)" +1 
2 " 2 2 "(/i— l)!(/i+1)! ’ "* U ‘ 

o 0 can be obtained directly from the initial expression (since P 0 = 1) 
and is equal to □ 

Problem 4.6 Expand (5—4x)~* as a series of Legendre polynomials. 

Solution. The straightforward approach using Theorem 11 will lead 
to evaluating integrals of the form j_' ,(5—4x) -i P n (x) dx and this does not 
seem to be a very simple task. It is, however, worth checking in the case 
of an expression of the form (a + bx)"* whether it is a particular case of 
the left-hand side of equation (4.17). In the present example 

(5—4x)‘* = (4 —4x+l)-* = Kl-2^x+i)-*. 

This suggests applying equation (4.17) with t = giving 

(5-4x)"* =| X 2-"P„(x). □ 

It = 0 


4.3 Fourier-Bessel Series In the present section we examine the possi¬ 
bility of expanding a function on [0.1] as a series of terms proportional to 
where 7 v (/ vn ) = 0, and the relevant expansion theorem is: 


Theorem 12 If fix) is piecewise smooth in [0,1] then the series 

00 

Z tfn-U/wr*), where 
n * 1 


a 


n 



xfix)Jfj yn x)dx. 


converges to /(x) at all points at which / is continuous and to 
i[/( x + 0)+/(x —0)] at points of discontinuity. The series converges 
to zero at x = 1 for all v and for v > 0 converges to zero at x = 0. The 
particular form for a n may be deduced from equations (4.3) and (4.26). 

Problem 4 J Obtain the representation in [0,1] of unity as a series of 
functions J 0 (j n x) where J 0 (j n ) = 0. 

Solution. This is a particular case of Theorem 12 with v = 0 and / = 1 
and thus 

a -~jk[ XJAx)dx - 
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The integral can be evaluated by applying equation (4.22) with v = 0 
and A = j n , giving 


a 


it 


2 

JJX) 


C V.ME 


= i _ 
= 0 — 


2 

JnW 


□ 


Problem 4.8 If p n are the roots of 7 0 (4p) = 0, obtain an expansion in a 
series of J 0 (p n x) of the function /(x) defined on [0.4] by f(x) = 1, 
0 < x < 2, /(x) = 0, 2 < x ^ 4, /(2) = f. 


Solution. Defining a new variable y by jx gives an expansion problem in 
[0,1] in the new variable. The j m of Theorem 12 have now to be replaced 
by 4/i n and thus the coefficient a n of J 0 (j 0n y)[ = J 0 (p n x)] is 


jm £ d ' - 7W> £ 

(equation (4.22), with v = 1, A = 4 p n ). 


□ 


The only cases where the coefficients can be evaluated are those when 
the integrals occurring are such that the equations (4.22) or (4.23) can be 
used. There are some other special cases where the coefficients take on a 
comparatively simple form, but these cases require a detailed knowledge 
of the properties of Bessel functions. Further expansion problems for 
Fourier-Bessel series and their generalisations occur in Problems 5-2, 
5-4,5-5,5-6, 5-9 and 5-10. 


4.4 Sturm-Liouville Theory It is effectively shown in Problem 2.1 that 
the set of functions sinwrx/c(n = 1,2,...) is generated by determining A 
such that X" = — kX has non-trivial solutions in [0, c] vanishing at the 
end points of the interval. This type of problem is known as an eigen 
value problem, the appropriate values of A are the eigen values and the 
corresponding solutions are known as the eigen functions. Another 
example is provided in Problem 2.3, where it is shown that the functions 
sin (2n-\)nx/2l are eigen functions for the above equation in [0,/] when 
the conditions at the end points are X(0) = X'(l) = 0; the corresponding 
eigen values being (2n — l) 2 jt 2 /4/ 2 . 

Both the above examples are special cases of a general class of eigen 
value problem, the Sturm-Liouville problem, which is such that an 
arbitrary function can be represented in terms of the solutions to the eigen 
value problem. 

The Sturm-Liouville problem in its most general form is the determina¬ 
tion of the eigen values and the corresponding eigen functions in [a,b] of 
the differential equation: 
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d_ 

dx 



+q<j)+kp<}) = 0, 


(4.28) 


where p, q, p are given functions and the boundary conditions at the end 
points are either 


dtp 

ix 1 (l>+a 2 — = 0, x = a. 


(4.29) 

dd> 

M + ^=0> x = b. 


(4.30) 

k^W+k^ib) = 0 = + m 2 (;j 

9... 

(4.31) 


The first class of boundary conditions are known as unmixed boundary 
conditions and the second class are normally called periodic boundary 
conditions. If p # 0 in [a,b] then the problem is said to be regular. 

The principal results of Sturm-Liouville theory may be summarized in 
the following theorem. 


Theorem 13 Let p. p', qp and (pp)" be continuous real functions of x on 
[a, b ] with p > 0, p > 0 and the constants a,, a 2 , /?,, P 2 , k v k 2 , m,, m 2 
be real. Then the Sturm-Liouville problem has an infinite number of eigen 
values A n all real and not more than a finite number being negative. The 
eigen functions (j> n , <f> m corresponding to different eigen values A n , k m 
are orthogonal with weight p. 

Any continuous function /(x), with a piecewise continuous first deriva¬ 
tive, can be represented in a < x < b as an infinite series of the form of 
equation (4.2) with the coefficients defined by equation (4.3). Furthermore 
if the function/ satisfies the same boundary conditions as those of the 
Sturm-Liouville problem then the series converges to f in [a. />]. 

For the case of unmixed boundary conditions (i.e. equations (4.29), 
(4.30)) there cannot be two linear independent eigen functions correspond¬ 
ing to the same eigen value. Also if the additional conditions a, a, < 0, 
P 2 P t > 0, are imposed then k n 2* 0. 

The orthogonality conditions also hold for the unmixed problem 
with p(a ) = 0 provided that the equation (4.29) is dropped and with 
p(b) = 0 provided that equation (4.30) is dropped. 

Note. The Legendre and Fourier-Bessel expansions of the previous 
sections are particular examples of irregular Sturm-Liouville problems 
and the orthogonality of these sets of functions is a direct consequence of 
the above statement concerning orthogonality for the particular case when 
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p(a) and /or p(b) are zero. The boundary conditions for the Sturm-Liou- 
ville problem for the Legendre expansion are that </> is bounded at x = 1, 
— 1 and those for the Fourier-Bessel expansion are that (j> is bounded at 
x = 0 and vanishes at x = 1. 

The set sin nnx/c clearly corresponds to the case p = 1 = p, q = 0 in 
[0,c] with x 2 = p 2 = 0 and the sets cos nnx/c. sin(2n — \)nx/c, 
cos(2n- l)7rx/2c, are also generated by the same equation with different 
sets of values of a, p [for cos nnx/c, a,=/*,= 0; for sin(2n- 1 )ttx/2c, 
a 2 = Pi = 0;forcos(2n—l)rrx/2c, a, = P 2 = 0]. The various expansions 
quoted in Chapter 1 are particular cases of Theorem 13. 

Problem 4.9 Find the Sturm-Liouville problem in [0,2c] whose eigen 
functions are cos nnx/c, sin nnx/c. 

Solution. These functions satisfy X" + ?.X = 0 with A = n 2 n 2 /c 2 , cos nnx/c, 
sin nnx/c do not satisfy the same mixed conditions at x = 0 and x = c 
but satisfy the periodic conditions <f>(0) = <j>(2 c), 4>'(0) = <j>'(2c), and we 
thus examine the Sturm-Liouville problem for the above equation with 
these conditions. For A # 0 the general solution is A cos A*x + B sin A^x 
and the boundary conditions give 

A = A cos 2cA 4 + B sin 2cA*, B - B cos 2cA* — A sin 2cA*. 

These equations will only have non-zero solutions if their determinant 
vanishes and this gives cos 2cA* = 1 and hence A = n 2 n 2 /c 2 . If this value 
of A is substituted into the equations they degenerate, and hence any values 
of A and B are possible and in particular cos nnx/c, sin nnx/c are both 
eigen functions, illustrating that the result concerning the non-multiplicity 
of eigen functions is not valid for periodic conditions. 

For A = 0 the general solution is Cx+D and the conditions give 
D = 2Cc + D,C = C, and hence C = 0 and in this case there is only one 
eigen function, i.e. cos nnx/c (n = 0). □ 

Problem 4.10 Determine the eigen functions and eigen values for 
d 2 <p/dx 2 +2<f> = 0, 0 < x < 1, with $(0) = 0 and h<f>(l)+<f>'(l) = 0. 

Solution. For A # 0 the general solution is A cos A*x -t- B sin A*x and 
as </»(0) = 0 we have that A = 0; the condition at x = 1 gives 
hB sin A 4 +A*# cos A* = 0. In this case the product p t P 2 of Theorem 13 
is equal to h and hence, for h > 0, A will be positive thus A* will be real. 
This can also be deduced directly by writing A = — u 2 which gives 
hB sinh u+Bu cosh u = 0, and, if B is not identically zero, tanh u = —u/h 
and for h > 0 this is not possible. If h is negative then, as tanh u is a 
monotonic increasing function of u, there may be one solution. The slope 
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of the tanh curve, however, is less than unity and, as both sides of the 
equation are equal for u = 0, it thus follows that the curves will not 
intersect for ft > -1. Thus for ft < -1 there will be one negative eigen 
value A = — u 2 , where u is the positive root of tanh u = — u/h and the 
corresponding eigen function is proportional to sinh ux. 

The remaining eigen values are given by A B = r 2 where the v n are the 
positive roots of ft tan v = —v. (It is obvious from sketching the graph of 
tant) that there will be an infinite number of solutions of this equation; 
this is also implied by Theorem 13.) The eigen functions are proportional 
to sin v n x. 

The general solution for A = 0 is <j> = Cx + D, and since <p vanishes at 
x = 0 it follows that D must vanish and the condition at x = 1 gives 
Cft + C = 0. Thus for the case ft = - 1, A = 0 will be an eigen value with 
the corresponding eigen function being proportional to x. D 

Problem 4.11 Determine the eigen values and eigen functions for 
x 2 d 2 (p/dx 2 + 3xd<(>/dx + ?.<!> = 0, in [l.e] with <p vanishing at the end 
points of the interval. 

Solution. This is not in Sturm-Liouville form but can be rewritten in 
that form with p = x 3 , p = x. The first step is the determination of the 
general solution and the form of the equation suggests considering <f> to 
be of the form x" and this gives (p + 1) 2 + A- 1 =0. Thus for A # 1 the 
general solution is 

<{> = y4x ,,-A)Vi-1 + Bx -<1-A,V4 _1 . 

The condition at x = 1 gives A = —B and application of the condition 
at x = e gives A sinh(l — A) 1 = 0. Thus (1 — A)* = inn and A = 1 +n 2 7t 2 , 
n = 1,2,.... and the eigen functions are proportional to x“ , [x" m -x"""']. 

The above analysis does not show that A = 1 is not an eigen value and 
this case has to be treated separately. For A = 1 it can be verified that 
independent solutions of the equation are x -1 and x" 'In x, and thus the 
general solution is x - l [C + D In x]. 

The condition at x = 1 gives C = 0, and from = 1 it follows that 
D = 0. Hence A = 1 is not an eigen value. □ 

Problem 4.12 Obtain, for the case ft > 0, an expansion of unity in [0,1] 
in terms of the eigen functions of Problem 4.10. 

Solution. The eigen functions are proportional to sin» B x where 
ft tan v =: —v and thus, by Theorem 13 and equation (4.3), 

n n * 

QC 

1 = X a n sin v n x, 

n = 1 
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where j‘ sin v n x dx = aj‘ sin\ x dx cos 2r„ x) dx, hence 

a n = 4( 1 - cos p.)/(2 - sin 2v n ). □ 


In using equation (4.3) we have assumed the result that the eigen 
functions of a Sturm-Liouville problem are orthogonal, and in this 
simple case we shall give a direct verification of this general result. 
Integrating by parts, we have 


sin v n x sin v m x dx = — sin v m cos v n + — T cos v 
Jo v n J 0 


x cos v_ x dx 


= — sin v cos v.—y sm v cos v 

\) m n yi n n i 


Hence 


v 2 p . 

J q sm t’ n xsin t^ 


X dx. 


( 1 V m\ f ' • • , COS » COS V, r v 1 

<) Jo sint, » xs,nt; ™ x</x =- f' - 2 [ tan ^-7tanp„J. 

As v „> v m both satisfy (tan v)/v = — l//i the right-hand side vanishes and 
hence, for v m # v n , so does the integral on the left-hand side. 


EXERCISES 

1. Expand x 5 in [ — 1,1] as a series of Legendre polynomials. 

2. Expand as a series of Legendre polynomials /(x) defined by: 

f(x) = 0, -1 < x < a, fix) = 1, a < x < 1. 

[Hint: Use (2n+ l)P„(x) = P n+ ,(x)-P;_ ,(x).] 

3. Obtain the Fourier-Bessel expansion of x 3 on [0,1]. 

4. Obtain the equation for the eigen values of d 2 y/dx 2 + ky = 0. subject 
to the conditions y(0)+/(0) = 0 = 2y(«)+3/(4 
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Chapter 5 

Solution of Boundary-Value Problems by Means 
of Generalized Fourier Expansions 

It is possible to use the results of Chapter 4 to extend the methods of 
Chapter 2 to problems when the functions a v a 2 , a 3 are not constant and 
when a linear combination of u and du/dx is prescribed on x = 0 and 
x = c. Problems where the equation has the form of equation (2.2) and 
a linear combination of u and du/dx is required to vanish at x = 0 and 
x = c will be termed homogeneous problems and, as in Chapter 2, the 
cases of homogeneous and non-homogeneous problems will be considered 
separately. 

5.1 Homogeneous Problems The technique is a simple generalization 
of that of §2.2 A product solution X Y is assumed and separation of 
variables used to determine one of the functions (which will be assumed 
to be X) as the solution of an eigen value problem and the principle of 
superposition is then used to obtain a general representation. If the eigen 
value problem is of a type such that any function can be expanded in 
terms of the eigen functions then the boundary conditions and equation 
for the corresponding Y can be obtained as for Fourier coefficients. It is 
possible to avoid some of the manipulative work of separation of variables 
in special cases by remembering the general form of product solutions of 
certain equations. It is shown in §4.1 that the general form of product 
solution of Laplace’s equation finite for all 6 in a spherical region is 

£ iA n r" + B n r~ n ~ l )P B (cos0). 

n - 0 

Another useful result which follows from the analysis of §4.1 is that 
[A J 0 iXp) + B V 0 (;.p)][C cosh kz + D sinh Az], where y o is the second solution 
of Bessel’s equation, satisfies the axisymmetric Laplace equation. 

Problem 5.1 Find a solution u of Laplace’s equation, valid outside the 
sphere of radius a, tending to zero at infinity and such that on the sphere 
u = cos 6, where 0 is the angle between the radius vector and a fixed 
direction. Obtain also the solution to the problem when the condition on 
the sphere is u = ia 2 + f 2 -laf cos0) -i , / > a. 

[The first problem corresponds to determining the total electrostatic 
potential when an earthed sphere is placed in a field of uniform intensity. 
The second problem arises when the uniform field is replaced by that of 
a point charge at a distance / from the centre of the sphere.] 
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Solution. If spherical polar coordinates (r, 0) are chosen with origin at 
the centre of the sphere, then the analysis of §4.1 shows that an appropriate 
expansion to choose for u is 

«= f (V"+cosfl). 

n = 0 

This form automatically satisfies Laplace’s equation and is finite for all 
real 0. In order to satisfy the condition at infinity A n = 0 and the condition 
on r — a gives 

00 

£ B a a~ M ~ l P(oos0) = cosO. 

n = 0 

Thus we need to expand cos0 as a series of Legendre polynomials, but 
cos0 = P,(cos0), hence B n = 0, n # 1, B t = a 2 . 

To solve the second problem we require the expansion of 

(a 2 + / 2 - 2a/ cos 0) " * = / " 1 [ 1 + ( fl / f) 2 - 2(a//)cos 0] " * 
and equation (4.17) shows that 

[ 1 + (a//) 2 - 2(a//)cos 0] "* = £ (a/J') a P n (cos 0), 

n = 0 

hence B n = a 2 " +1 //" + '. The solution is thus 

u = ~ £ j(j r ) p n(* os0) = j [r 2 +o 4 // 2 -2(a 2 r//)cos0]-*, 

n = 0 

from equation (4.17). The quantity in square brackets is the inverse 
distance from the point (r,0) to (a 2 //,0) and the above form of solution 
is that obtained directly in the electrostatic problem by the method of 
images. q 

Problem 5.2 Solve * 

d 2 u 1 8u d 2 u 

= 0$r<a,t>0, 

subject to the conditions 

u(0,r) bounded, u(a, t) = 0, f > 0, u(r,0) = (1 -r 2 /a 2 \ du(r, t)/dt = 0, 
t = 0,0 ^ r < a. 

(This corresponds to determining the displacement at any time of a 
circular membrane of radius a which is fixed round its edge, the initial 
displacement is 1 — r 2 /a 2 and the membrane is initially at rest.) 

Solution. We assume a solution of the form F(r)G(r) and the method of 
separation of variables gives 
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dr 2 


r dr 


0 , 


d 2 T . 

-rJ-+;. 2 T = 0. 

dt 2 


The equation for F is Bessel’s equation of order zero, in order that u(0, t) 
be bounded F must be J 0 (/.r) and the condition on r = a gives J 0 (/m) = 0. 
Thus Xa = j 0n and to avoid unduly complicated expressions we drop the 
suffix zero. Thus the appropriate form of solution is 

00 

h = Z ^ 0 (fn cos/, t/a + B n smj n t/a]- 

n = 0 

Since du/dt = 0 for t = 0 it follows that B n = 0 and thus A n is the coeffi¬ 
cient of J 0 (j n r!a) in the expansion of 1 —r 2 /a 2 in terms of J () (j n r/a) or, 
equivalently, of J^{j n x) in the Fourier-Bessel series in [0,1] for (1 — x 2 ) 
(this follows on writing r/a as x). Theorem 12 gives 

A„ = -j^jr ) | (l-x 2 )x J 0 {j n x)dx. 

The first part of this integral has already been evaluated in Problem 4.7 
and the second part can be evaluated as follows: 


| dx=j | x 2 ^[x y t (/„ x)] dx. 


IMA 

K in 


f. 


X 2 J,UnX) ^ 


from (4.22), with v = 1, 


m_2 

in fn 


on integration by parts, 
j 0 dx, 

from (4.22), with v = 2, 


in in 

Thus from the above result and that of Problem 4.7, 



this can be further simplified by using equation (4.24) and the fact that 
J 0 (j n ) = 0 and the final result is A n = 8/j 2 /,(/,). □ 


Problem 5.3 Solve du/dt = d 2 u/dx 2 , 0 < x < 1, t > 0, under the 
conditions u(0,f) = 0, du/dx+hu = 0, x = 1, t > 0; u(x,0) = 0, 

0 < x < 1, h being a positive constant. 

Solution. The occurrence of d 2 u/dx 2 suggests a Fourier series in x but 
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the condition on x = 1 is not appropriate for Fourier series. We therefore 
assume product solutions of the form X(x)T(t) and carrying out the normal 
separation of variables procedure gives 

X"+kX = 0, T' + kT = 0. 

The conditions on x=0 and x = 1 give JT(0) = 0, X'(l)+hX(l) = 0 
and thus the eigen value problem for X is precisely that of Problem 4.10. 
Hence, using the notation of that problem, 

00 

« = Z A n exp(—^ t)sin v n x, 

n = 1 

and on setting r = 0 it follows that the A n are the generalized Fourier 
coefficients of x. The weight p in this case is unity and thus 

A n P sin 2 v n xdx = P xsinuxdx, 

i.e. v n (2v n - sin 2 v n )A n = 4(sin v n - v n cos v n ). □ 

Problem 54 Solve 


d 2 u 1 du d 2 u 
d? + rfr + d? = °’ 


0<z<h,0<r<a, 


subject to the conditions u(0, z) finite, u(r,0) = 1, du(r,z)/dr+yu = 0, 
r = a; du(r, z)jdz + fiu = 0, z = h, where /?, y are positive constants. 

[One physical interpretation of this problem is the determination of 
the steady state temperature in a cylinder when one of its plane faces is 
maintained at unit temperature while a radiation condition is imposed 
on the other surfaces.] 

Solution. The application of separation of variables to this equation 
has been considered in §4.1 and appropriate product solutions satisfying 
the finiteness condition are 7 0 (Ar)[/l cosh A(z-/i) + B sinh A(z-/i)]. Since 
homogeneous conditions are imposed on z = h the algebra is simplified 
slightly by using hyperbolic functions of z-h. The boundary condition 
on r = a gives 

j r + 7 7 0 (Ar) = 0, r = a, 
i.e. XJ' 0 ()M)+yJ 0 (Xa) = 0, 

where the prime denotes derivative with respect to the argument. 
[A'o(z) = — J x (z)].Thus ka will have to satisfy the equation^ J' 0 (pi) + 6J 0 (jj) = 0 
where S = ya. It can be shown that there are an infinite number of solutions, 
all positive, of this equation and of its generalization n J^(p) + d = 0. 
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It can also be shown that any function /(x) continuous in [0,1] can be 

OO 

represented by a series of the form £ a n J s in vn x) where // vn are the roots 

n - 1 

of the above equation. The coefficients a n are given by (G. N. Watson Theory 
of Bessel Functions , C.U.P., 1952, Chapter 18) 

2Pvn I* t/(t)7 v (^,t) dt 

a ” 

This more general type of series is known as a Dini series and is an obvious 
generalization of the Fourier-Bessel series (which corresponds to the case 
5 -* oo.). For the particular case v + 5 = 0 the correct expansion for a 
function as a Dini series is (G. N. Watson, Theory of Bessel Functions, 
C.U.P., 1952, Chapter 18) 

Z a n JM,x)+2(v+l)x v f ' t v+l f(t)dt. 

n = 1 

Reverting to the original boundary value problem it is now seen that the 
appropriate representation is 

00 

Z y (M, r/a)\_A n cosh \xfz - h)/a + B n sinh n n (z - h)/a ], 

n = 1 

where the suffix zero is dropped to avoid unnecessarily complicated forms. 
In order to satisfy the condition at z = h, fi n B n + fiaA n = 0 and the 
condition on z = 0 then gives 


Z Jo(M H r / a )A K [p°sh n n h/a+(fia/n ll )smh n n h/a] = 1, 

n = 1 

and A n can now be calculated from the quoted result for Dini series, the 
integral being evaluated as in Problem 4.7. The final result is 

A„[n„ cosh n n h/a + Pa sinh n n h/a] = 27, (n n )/J 2 0 iF „)+•/?(/*„)• □ 


Problem 5.5 


Solve 

d 2 u 1 du du 


dr 2 + r dr dt’ 


C 

t > 0, 0 < r < 1, 


under the conditions u(0,t) finite, du(r,t)/dr = 0, r = 1, t > 0; u(r, 0) = r, 
0 < r < 1. 


Solution. Following the standard separation of variables procedure 
shows that there are separable solutions of the form 
{ A 7 0 (Ar) + B y o (Ar)}exp -A 2 f. 

The finiteness condition gives B = 0 and the condition on r = 1 gives 
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7g(/l) = 0. Thus the set of values for A are P n where P n are the non-zero 
roots of Jg(x) = 0 and A can also be zero. The general form for u is 

*o+ t a n JJP 9 r)pxp(~P 2 n t) 

n * 1 

and the condition on t = 0 gives 

a 0+ Z a n J 0 (p K r) = r. 

It = 1 


The expansion problem is thus a particular case of the Dini series of the 
previous problem and it is also one of the special cases when v + S = 0. 
Thus 

In this case the integral for a n cannot be simplified. 



Problem 5.6 Solve 


c ( du\ d 2 u 
dx\ dxJ dt 2 ' 


0 < x < l, t > 0, 


under the conditions u(0, f) bounded, u(/,r) = 0, t > 0; u(x, 0) = /(x), 
where/isa given continuous function ofx, du(x, t)/dt = 0,r = 0,0 < x < /, 
(This problem arises in the investigation of the small oscillations in the 
vertical plane of a heavy chain displaced from rest). 


Solution. This equation is not one for which the product forms of 
solution can be instantly recognized, but writing u as X(x)T(t) gives 

dX\ d 2 T 

dx) +XX = °’ 1S +XT =°- 

X is required to be finite at x = 0 and zero at x = / and the problem of 
determining X is a singular Sturm-Liouville problem. A change to the 
independent variable y(= x*) gives 




+ 4JLX 


= 0 , 


This is Bessel’s equation of order zero in the variable 2and, reverting 
back to x, the solution finite at x = 0 is J 0 ( 22*x*). The condition at x = / 
gives 2(A/)* = j n where the j n are defined as in Problem 5.2, and thus the 
appropriate form of solution is 


“ = £ Ucos^+Vin^kc ijtxm. 

It- 1 v J 
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The condition on du/dt at t = 0 shows that B n = 0 and the other condition 
at t = 0 gives 

/(*)= t \JoWWl 

It - 1 

This equation is not in a form such that the Fourier-Bessel theorem can 
be applied immediately but this can be remedied by writing ( x/f) i as z 
and it is then seen that the A n are the coefficients in the Fourier-Bessel 
expansion of f(z 2 l), 0 < z < 1, i.e. 

j ° u - i)dz - 4I>^) *■ n 

Problem 5.7 Solve (l 2 — x 2 )d 2 u/dx 2 — Ixdu/dx — d 2 u/8t 2 , 0 < x < l, 
t > 0, under the conditions u(0, t) = 0, u(l, t) bounded, t > 0;u(x,0) = x 3 , 
du(x,t)/dt = 0, t = 0, 0 < x < /. (This problem occurs when investigat¬ 
ing the small transverse vibrations of a string rotating about a vertical 
axis.) 

Solution. There will exist product solutions of the form X(x) T(t) where 
(l 2 - x 2 )X" - 2xX’ + AX = 0, T"+XT = 0. 

The equation for X is Legendre’s equation in the variable x/l but only 
positive values of the argument are considered. If X* is defined by Af(x) 
for x ^ 0 and by - Af(-x) for x ^ 0 then X* satisfies the same equation 
as X, is continuous for all x with |x| ^ I. X* is also bounded at x = +/ 
and thus must be a Legendre polynomial of the variable x/l. Further, 
since X* is odd in x it must be a Legendre polynomial of odd order, i.e. 

for x > 0 X* = X = P 2n _, ( x/l) and A = 2n(2n-l). n = 1. The 

general form of solution is thus 

u = I M 21 .- 1 cos [2«(2n — l)]*t + B z j sin[2n(2n — l)]*f }P 2 „- ,(jc//), 

the condition on du/dt gives B 2n _ { = 0, and from the value of u(x,0) 
it is found that 

* 3 = I Au-iPu-iWl 

It - 1 

^ 2 n -1 cou ld be found by an extension of the analysis of Problem 4.4 
but it is easier in this case to use the result 2P 3 (x) = 5x 3 - 3x, i.e. 
x 3 = ±[2P 3 (x)+3P,(x)] and 

x 3 = i/ 3 [2P 3 (x/0+3P 1 (x/0], 

i.e./I, = 3/ 3 /5,/4 3 = 2l 3 /5,A 2n _ 1 =0, n # 1,2. 


F 


□ 
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Problem 5.8 Obtain a solution of Laplace’s equation within the sphere 
r ^ a, finite everywhere within this region, equal to unity on the sphere 
for 0 < 9 < a and zero on the remainder of the sphere (r, 0 being the 
usual spherical polar coordinates). 


Solution. The appropriate form of solution to use is that of Problem 4.1, 

in this case however the finiteness condition within the sphere requires 

00 

that B n = 0 and the appropriate form of solution is £ /t„r"P n (cos 0) 
with " = 0 


~ . __ . .. 1, 0 < 6 < a, 

Y A a P (cos 9) = . 

„ = „ 0, a < 0 < 7i. 

The A n can now be obtained from Theorem 11 on writing x as cos 6; 
this gives A n tf = \(2n + 1) j‘ P n (x) dx. For n = 0 the integral is equal 
to 1 — cos a and for n ^ 0 it can be evaluated by use of the recurrence 
relation 

(2n+l)P„(x) = F + 1 (x)-P'_,(x), 

where the dash denotes the derivative with respect to the argument. Hence 

A n a n = 2[ /> n+i( cosa )- p »-i( cosa )]’ n >1. □ 


Problem 5.9 Solve 

d 2 u 1 8u 6u 

-a? + -rfr - 7,’ 0 <z<h,a<r<b, 

under the conditions u(a, t) = u(b,t) = 0 ,r > 0;andu(r,0) = 1 ,a < r <b. 


Solution. Separation of variables shows that there will be product 
solutions of the form {A J 0 (Xr) + BY 0 (Xr)}exp{- X 2 t) and the conditions on 
r = a and r = b give 

A J 0 (Xa) + B Y 0 (Xa) = 0, A J 0 (Xb) + B Y 0 (Xb) = 0. 

For non-zero A and B the determinant must vanish and thus 
J 0 (Xa)Y 0 (Xb)-J 0 (Xb)Y 0 (Xa) = 0, 


which is the equation satisfied by A. 

An equivalent method of calculating A would be to use the fact that 
there will be product solutions of the form R(r)exp( — X 2 t) where 
R" + R'/r+X 2 R = 0 and R(a) = R(b) = 0. The problem looked at from 
this point of view is of standard Sturm-Liouville type and from this 
general theory it is deduced that there will be an infinite number of values 
for A (which we shall denote by <5 W , n = 1,..., oo) and the eigen functions 
will be orthogonal to each other with weight r. The nth eigen function is 
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The general form for u is thus 

X “n exp( -<5„ 2 0[y o (^ n r)y 0 (<5„ a)- Y 0 (S n r)J 0 (S n a)], 

n = 1 

x 

where the a n are determined by £ a n (j> n = 1. On using the orthogonality 

n — 1 

property of this gives 

J/4>n dr = j%<£ n (r) dr. 

The evaluation of the integrals in the above expression is somewhat 
complicated, but it can be completed by using various relations involving 
Bessel functions [equations (1) and (10) of G. N. Watson, Theory of Bessel 
Functions , C.U.P., 1952, Chapter 5]. The final results are somewhat 
complicated and hence will not be quoted. □ 

5.2 Inhomogeneous Problems The extension of the lethods of the 
previous section to inhomogeneous problems will now be considered. The 
first step is to determine the set of functions (j> n which would be most 
appropriate to use, these can be found from considering the corresponding 
homogeneous problem as illustrated in the previous problems. The 

X) 

form of the solution is thus £ Y n (y)(f) n (x) and the next step is determining 

n = 1 

the equation for Y n . The analysis will be restricted to equations of the form 

mhk M I) +0 » w " +fr>(y) lp + b ‘ (y) ^ + ‘w* - < 1 

on using the method of separation of variables it can be verified that the 
problem of solving the homogeneous problem reduces to a Sturm- 
Liouville problem and that p(x) is the weight function. The above equation 
is then multiplied by p<j> n and integrated with respect to x from x = 0 to 
x = c. Integration by parts twice on the left hand side will give that this 
side of the equation will involve derivatives with respect to y of 
Jo P u [* e - T(y) p<t>ldx~\dx and the prescribed boundary values. The 
right-hand side will be the known function J* p<j> n f (x, y) dx; thus a differen¬ 
tial equation will be obtained for the y r For the case when the boundary 
conditions are homogeneous the equation for the Y n can be obtained more 
simply by expanding/ in terms of the appropriate <f) n . Substitution of this 
expansion and the corresponding expansion for u into the differential 
equation and equating coefficients of (j> n will then give the differential 
equation for y n . 
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In elliptic problems, when conditions are prescribed for two values of 
x and two values of y and one of the two sets of conditions is homo¬ 
geneous then it is best to use an expansion which satisfies these conditions 
automatically. 

It should again be remembered that the simplest method, if possible, is 
to reduce the inhomogeneous problem to a homogeneous one. 


Problem 5.10 Solve 

d 2 u 1 du d 2 u 
dr 2 + r dr dt 2 


—<7 sin cot. 


0 $ r < a, t > 0 , 


under the conditions u = du/dt = 0 at r = 0 , u = 0 , r = a; u finite for 
all r, 0 s* r < a, aw not being a zero of J 0 (z). 


Solution. The boundary conditions are homogeneous and it follows 
(on considering the case q = 0) by direct comparison with Problem 5.2 

oo 

that the appropriate expansion is u = £ ^oOi/AO^MO where we use 

n = 1 

the notation of Problem 5.2. The appropriate expansion for the right- 
hand side of the equation is obtained from Problem 4.7. Substitution in 
the equation and equating coefficients of J 0 (j„r/a ) gives 




2 q sin wt 
Jn J M ' 


The conditions on u, du/dt at t = 0 show that i/^ n = ^ = 0, t = 0, and, 
provided j\ ^ co 2 a 2 , the solution of the above equation satisfying these 
conditions is 

= 2 qa 2 [/„sin wt-jo sint] Q 


Problem 5.11 Solve d 2 u/dx 2 = du/dt, 0 < x < a, t > 0, under the 
conditions u(x, 0 ) = 1 , 0 < x < a; u{a,t) = 1 , du/dx = hu, x = 0 , t > 0; 
h>0. 

Solution. The simplest approach is to look for a solution which satisfies 
the conditions at x = 0 and x = a. Such a solution is (1 + hx)/(l +ha) 
and writing u as (l + /ix)/(l + ha)+v shows that v satisfies the same equa¬ 
tion as u and v satisfies homogeneous conditions on x = 0 and x — a\ 
this problem can then be solved by the previous methods. 

We shall, however, consider the application of the direct technique 
described above. The homogeneous problem is that when the boundary 
condition on x = a is u = 0, and by comparison with Problem 5.3 the 
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appropriate expansion for the homogeneous problem is 
Y, 2f ll (x)exp(— X n t) where X" + X H X = 0, and X' = hX , x = 0; X = 0, 

n = 1 

x = a. 

The solution for X n satisfying the condition on x = a is £ n sin a*(x — a) 
and the condition on x = a gives A*cos = — /7 sin ft a. Thus X*a = w n , 

where vv n are the roots of Ziatanvv = — w. The form of solution to be 

00 

sought is thus u = £ i/^(f)sin w w (x — a)/a and 

n= 1 

il/ n sin 2 w n {x — a)/a dx = u(x, t)sin w n (x — a)/a dx . 

It can be shown by integration by parts twice and using the boundary 
condition on x = 0 , a and the definition of that 

( a d 2 u . *v n (x —a) w n w 2 f a . w„(x-a) 

J 0 dx 2 a a a 2 Jo a 

Hence multiplying both sides of the partial differential equation by 
sinvv w (x-a)/a and integrating from x = 0 to x = a gives, on using the 
definition of \// n in terms of u. 


dil/ w 2 / 

jf- = - * - w nj a J o S| n 2 w(x - a)/a dx 


Thus 


= —S-i b - 
a 2 Wn 


4w 2 


" o 2 ( 2 w n —sin 2 w n )’ 


-sin 2 w n ) 

00 

The condition on r = 0 gives 1 = £ ^ n (0)sin w n (x - a)/a, i.e. 


r. 


n = 1 

. W 0 ) 


f . 

sin 2 — 2 - dx 

Jo a 


or 


Hence 


a at// n ( 0 K 2 w -sin 2 w ) 

—(cosw„-l) = - " — 2 - 

4w. 


w 


A = 


4 cos w 


( 2 w n —sin 2 w n ) 


Problem 5.12 Solve 


d 2 u 1 du d 2 u 

3 7 + -rTr + J? ’ 0< r <‘' °<’<*. 
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under the conditions 

u(0, z) finite, u(a,z) = u v 0 < z < h\ u(r, 0) = u 2 , u{r<h) = u 3 where 
u x ,u 2 ,u 3 are constants. 

Solution. The boundary conditions on z = 0 and z = h are simple so 
that it is worth considering the possibility of obtaining homogeneous 
conditions at z = 0 and z = h. The function u 2 (l — z//i) + u 3 z/h satisfies 
the conditions at z = Oandz = K and ifn is written as u 2 ( 1 -z//i)+u 3 z//i+ v 
then v will satisfy the same equation as u; v vanishes at z = 0 and z = h 
and is equal to u x — u 3 z/h — u 2 (\ - z/h) on r = a. 

This problem for v is of the standard type encountered in the previous 
chapter and can be solved by assuming a Fourier sine series representation. 

An alternative method of solution would be to attempt to obtain a 
series representation for u which would automatically satisfy the condition 
on r = a. One such solution is u x ; writing u as w + t^ shows that w satisfies 
the same equation as u, vanishes on r = a and is equal to u 2 — u , and 
u 3 — u x on z = 0 and z = h respectively. The problem for w can be 
solved by use of separation of variables and by analogy with previous 
examples it can be found that the appropriate expansion is the Fourier- 
Bessel series. □ 

In examples of this type where two approaches are possible, it is pre¬ 
ferable to use the one with least complications. In particular, if one set of 
conditions are homogeneous one should use the approach which auto¬ 
matically satisfies these conditions. Thus if u 2 and u 3 had been zero the 
sine series approach would be the less complicated, and for u { = 0 use 
of the Fourier-Bessel method would have advantages. There is the 
additional point that integrals involving Bessel functions are usually not 
as simple as those involving trigonometric functions, and this could mean 
that even for non-homogeneous conditions it might be better to use 
trigonometric series where possible. 


EXERCISES 


1 . 


Solve 


d 2 u 1 du d 2 u 
d? + 'i : dr + d? = °’ 


0 ^ r < a, z > 0, 


subject to the conditions w(0,z) finite, u(a,z) = 0, z > 0, t<(r,0) = 1, 
0 ^ r < a, u->0asz->oo. 

2. Obtain a solution of Laplace’s equation finite everywhere within the 
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sphere r = a and equal to cos 3 0 on the boundary of the sphere, r,0 being 
the usual polar coordinates. 

Obtain also the solution of Laplace’s equation valid outside the sphere, 
tending to zero at infinity and taking the same values on the boundary. 


3. Solve 


d 2 u 1 du d 2 u 


subject to the conditions du/dr = 0, r = a, w(0, r) finite, t > 0;w(r,0) = r, 
du/dt = r 3 , t = 0,0 ^ r < a. 


4. Solve du/dt = d 2 u/dx 2 , 0 < x < Z, t > 0, subject to the conditions 
du/dx = /iu, x = 0, du/dx = — hu, x = /, t > 0; u(x,0) = 1, 0 < x < /, 
h being a constant. 

5. Solve the boundary value problem of the previous exercise when the 
conditions on x = 0, x = / are du/dx = hu — hu x , du/dx = —hu + hu l 
respectively where u, is a constant. 
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Answers to Exercises 


Chapter 1 

1. (iv) n 2 /S, 7t 4 /96. 

( — lf + l sin mx 


2. x-2 l 


, 0 ^ x < 7c; 


cos 


A _ 2 4 y ( —iy +1 cosnx 

2 * n n L 1 - 


4m 2 


, 0 < x s* c. 


4. 


Z 


[(— 1 ) n +1 (2n — 1 )n exp(c) —2c] cos(2n — 1 )nx 


[2(— l) n+1 cexp c4-(2n — l)7t] sin(2n — l)7rx 


n = 1 
oc 

<Z 

n = 1 

2 4 A cos2mx 


(2m — l)V + 4c 


2c 


[(2n-l)V + 4c 2 ] 


2c 


5.- y 

% n L. 


n = 1 


4n 2 -l 


Chapter 2 

4pa y sinh[(2n-l)7ry/a]sech[(2n-l)7rfc/a]sin[(2n-l)7rx/a] 
n 2 (2n — l) 2 

-i • /i 1 • n 1 2 V P 2n ° os 2n0 / _ 

2. p sm 0, -p sin 0 4 - ) —— = —— (see Problem 1.6). 

2 n n L 4rr — 1 

n = 1 

4 sinh[(2n-l)jt>»/2a]cosech[(2n-l)7rfe/2a]sin[(2n-l)7rx/2a] 

(2n-1) ‘ 


3. 


A w 

;Z 


n ~ 1 


4. m — 1 is a solution of exercise 1. 
x 4 a 


5. a + (b-a) j-— Y sin[(2n— l)7tx//]exp[ — (2n— l) 2 7c 2 t// 2 ](2n — l)' 1 

n = 1 

21 °° 

-( b—a ) y ( — l)" + 1 sin(Mrtx/Oexp[ — nVf// 2 ]n _1 . 

71 Li 

6. £ ( -ir 

71 L-t 


7 - 2 Z 

»t = 1 


1+1 cos[(2m — l)irt/2/]sin[(2n — l)7tx/2/] 

[M7tsin t— /sin M7tt//] sin mttx 
(mV-/ 2 ) / ‘ 
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Chapter 3 

2. 2/J//? 4 +4. 

3. 2(cos x — cos 2x)/nx. 


2_r sin/?x r 

• *J r+vL 


/S 3 sin cot — fico cos cot 4- pco exp( 


- 0 * 0 ]#. 


(The first two terms can be evaluated explicitly using Tables of Fourier 
transforms (Erdelyi, A., ef a/. Tables of Integral Transforms, McGraw- 
Hill, 1954, vol. 1, p. 66, nos. (23) and (24)) and are equal to 
exp[ - x(ay2)*]sin[ft>t—x(co/2)*].) 


, -T 

7t J n 


cos Bx 

[exP( — 0 2 O—exp( — t)] dp. 


Chapter 4 

8 4 3 

1 _ P +-P +-P 

63 5 + 9 3 + 7^ 1- 

2. See Problem 5.8. 

3 - 2 y U 2 (j n )],using notation of Problem 5.2. 

nm\ J » J l Vi* 

4. (24-32)sin2*7t4-2*cos2*7r = 0. 

Chapter 5 

1. (The notation of Problem 5.2 is used in 1 and 3.) 

O V J o(j n r/a) _, . , , 

2 Z 77777 exp, - J - z/a) - 

2. f r 3 P 3 (cos 6 )+f r cos 9, fr~ 4 P 3 (cos0)+§r 2 cos0. 

3. 2a f J 0 (j n r/a)[j n 2 J fd m )co$j n tla+a*[j n Jfd n ) 

-2Jfj n )-]smd n tla-]ldlJ]{d n ). 

Y [A n sin X n l+h-h cos X n / ][cos X n x + h A~ 1 sin A n x] exp( - A 2 1 ) 

(X 2 n +h 2 )+2h ’ 

where A n are the roots of (A 2 -h 2 )tan A/ = 2/iA. 

5. Writing u as m 1 4- v gives that v is (1 — «,) times the solution of exercise 4. 
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Appendix 1. Some Simple Fourier Series 
sin(wrx/c) n{l — x) 


Z smynn: 
n 


2 c 


0 < x ^ c, 


*z 

it — i 

00 

3- Z 

n = ] 

x 

<■ Z 


cos(«7tx/c) _ — ln(2 sin nx) 
n 2c 


, 0 < x ^ c, 


( —1)" 1 sin nnx/c nx 

-= —, 0 $ x < c, 

n 2c 


(—1)" 'cos nnx/c 


-ln(2cosg). 

5 | 


0 $ x < c, 


X ^ c. 


s - z 


(—1)" 1 cos nnx/c 


~ n ih-£} 0 

7 - Z ^ - <£-£♦£} ° 

u = 1 x ' 


^ x ^ c, 


$ X $ c. 


O V (— 1)" - 1 sin nnx/c n 3 xf t x 2 \ n 

8 - Z—?-"M 1 -?/ 0 

If = l ' / 


s* x £ c. 
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Appendix 2. Table of Fourier Integrals 



fix) 

F(/l) = j exp( - ipx)f(x) dx 

J - <X> 

1. 

rye**-*) 

FUP) 

2. 

F(x) 

2nf(-fi) 

3. 

/(Axjexp ibx 

X- l F[(p-b)/X] 

4. 

(a 2 + x 2 )~ 1 

(n/a)exp—a\p\, a > 0 

5. 

exp— x 2 

7t*exp—j/? 2 


f(x) 

F,(/?) = J“ /(x)sin Px dx 

6. 

F,(x) 

W(P) 

7. 

r'f(xr l ) 

FMP) 

8. 

x” 'sinh ax cosech yx 

tan” '[tanf^roxy - ')tanh(^7iy “'/?)] 

9. 

x~ 1 exp—ax, a > 0 

tan ~ i a~ 1 P 

10. 

exp —ax 

P(« 2 +P 2 )~ l 


f(x) 

F C (P) = J* f(x)cospx dx 

11. 

F c (x) 


12. 

exp-ax 2 

^jr*a - *exp— p 2 /4ot 

13. 

cosh ax sech yx 

_ jcos(“7ty~ 1 a)cosh^7ty“ , /J 
cos Jiay -1 + cosh ny~ x P 
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